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ABSTRACT 

We present the results of a systematic study of mid-IR spectra of Galactic 
regions, Magellanic Hll regions, and galaxies of various types (dwarf, spiral, 
starburst), observed by the satellites ISO and Spitzer. We study the relative 
variations of the 6.2, 7.7, 8.6 and 11.3 /im features inside spatially resolved objects 
(such as M 82, M 51, 30 Doradus, M 17 and the Orion bar), as well as among 
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90 integrated spectra of 50 objects. Our main results are that the 6.2, 7.7 and 
8.6 /xm bands are essentiaUy tied together, while the ratios between these bands 
and the 11.3 /xm band varies by one order of magnitude. This implies that the 
properties of the PAHs are remarkably universal throughout our sample, and that 
the relative variations of the band ratios are mainly controled by the fraction of 
ionized PAHs. In particular, we show that we can rule out both the modification 
of the PAH size distribution, and the mid-infrared extinction, as an explanation of 
these variations. Using a few well-studied Galactic regions (including the spectral 
image of the Orion bar), we give an empirical relation between the I6.2/I11.3 ratio 
and the ionization/recombination ratio G^^jn^ x ^Tgas, therefore providing a 
useful quantitative diagnostic tool of the physical conditions in the regions where 
the PAH emission originates. Finally, we discuss the physical interpretation of 
the I6.2/I11.3 ratio, on galactic size scales. 

Subject headings: dust - HIT regions - ISM: structure - galaxies: dwarf, starburst 
- infrared: general 



INTRODUCTION 



The reprocessing of stellar light by dust in the infrared (IR) is widely used to probe em- 
bedded star formation. However, in the absence of other constraints, the physical properties 
which can usually be derived from an almost-featureless grain continuum emission is limited 
to global quantities, such as the dust mass and its average temperature. The ubiquity of 
numerous mid-IR aromatic features, in a wide variety of astrophysical objects and environ- 
ments, potentially provides more articulate diagnostics of the phys ical conditions. Indeed 



these features domi nate the mid-IR sp ectra of evolved stars fe.g. iBlommaer 



200 



3), as well as whole galaxies (e.g. Verma et al.l 2005 : 



: et al 



Kraemer et al.ll2006l ). the cool ISM fe. g. lAbergel et al.ll2005l:lFlagev et al.ll2006[|Povich et al 

Smith et al.ll2007l ). that have been 



2005: 



extensively observed by the Infrared Space Observatory (ISO)^ and are currently investi- 
gated with a higher sensitivity by the Spitzer Space Telescope. In our Galaxy, one third of 
the stellar light is reprocessed by dust, while this fraction can go up to 99% and higher, in 
starburst galaxies. At solar metallicity, roughly 15% of the cooling is radiated through the 
most powerful mid-IR bands, centered at 3.3, 6.2, 7.7, 11.3 and 12.7 /xm. 

Historically, these emission features were attributed to very small grains 10 A), tran- 
siently heated by single photon absorption, in order to account for the independence of the 
colo r temperature with the distance from the illuminating star, in several reflection nebu- 
lae (ISellgrenlll984l ). In parallel, the central wavelengths of these bands were recognized to 
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coincide with the vibrational modes of aromatic material (IDuley fc Williamslll98ll ). These 
features are now commonly attributed to the molecular modes of Polycyclic Ar omatic Hy- 
drocarbons (hereafter PAHs: lLeger fc Pugetlll984l : lAUamandola et al.lll985l . Il989l ). which are 
planar molecules made of 10 to 1000 carbon atoms, excited primarily by ultraviolet (UV) 
photons. With sihcate and carbon grains, PAHs are a main component of dust models 



flDesert et all Il990l : IPwek et all Il997l : IPraine fc Lil l200ll : IZubko et all 1200J) . Their absorp 



and quantum theory (in particular 


Desert et al. 


1990 


• 


Joblin et al. 


1992: 


Verstraete et al. 


2001 


: Li & Draine 


2001; 


Mattioda et al. 


2005a 




Draine & Li 


2007 




Malloci et al. 


200?!). In 



addition to being major radiative coolants of the interstellar medium (ISM), PAHs are re- 
sponsible for most of the photoelectric heating of the gas in photodissociation regions (here- 
after PDRs)_jJid__the_^n^^ their high cumulative surface area 
(e.g. iBakes fc TielensI Il994l : iHoUenbach fc TielensI 119971 ) . For the same reason, they pr oba- 



bly play an important role in grain surface chemistry (e.g. iTielens fc AUaraando 



2004 



a 



19871 ). In 



with solar 



our Galaxy, they contain 15 — 20% of the depleted carbon (jZubko et al.l 
abundance constraints). As a consequence, they are believed to be p art of the interstellar 
carbon condensation chain (ICherchneff et al.ll2000l : iDartois et al.l 120051 ). 



From an extragalactic p oint of view, the luni inosity of the 6.2 /xm feature can be used 
as a tracer of star formation (jPeeters et al.ll2004bl ). However, this tracer is biased by global 
parameters s uch as the ISM metallicity . Indeed, PAHs are under abundant in low-metallicity 
galaxies (e.g. iGalliano et al.ll2003l . [20051, l2008l : iDraine et al.ll2007l ) . There is a general correla- 



tion between the PAH-to-continuum intensity ratio and the ISM metallicity (IMadden et al. 



20061 : IWu et al.l l2006l : lO'Halloran et al.l l2006l : ISmith et al.l l2007h. and consequently be tween 
the IRAC8/im/MIPS24^m broadband ratio and the metallicity (lEngelbracht et al.ll2005l ). The 
ori gin of this trend has been attributed to radiative a nd mechanical destruction mechanisms 
bv lMadden et al.l fl2006l ) and lo^Halloran et al.l fl2006l ) respectively. Conversel y, from the de 



tailed modeling of the spectral energy distribution (SED) of nearby galaxies, IGalliano et al 



( 120081 ) showed that the PAH-to-gas mass ratio at different metallicities coincides with the 
relative amount of carbon dust condensed in the envelopes of low-mass stars, during the 
Asymptotic Giant Branch phase (AGB). This study suggests that PAHs are injected into 
the ISM by their progenitors, the AGB stars, several hundreds of Myr after the beginning of 
the star formation, when the gas has already been enriched by more massive stars. This delay 
corresponds to the time needed for AGB stars to evolve off the main sequence. Therefore, the 
delayed injection of AGB-condensed carbon dust into the ISM offers a natural explanation 
for the paucity of PAHs in low-metallicity environments. 



From a cosmological point of view, the large luminosity in these IR emission features 
coupled with the high sensitivity of Spitzer has allowed the detection of PAHs in distant 
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luminous infrared galaxies out to redshift z c::^ 2 (e.g. lElbaz et al.l l2005l : lYan et al.l l2005l : 
Houck et al.l l2005l ). Hence, understanding what controls the properties of the aromatic 
features on large scales is required, in order to properly interpret broadband surveys. 

The detailed characteristics of the mid-IR features, such as their shape, their cen- 
tral wavelength or t he intensity ratio between the different bands are known to vary (see 
Peeters et al.ll2004al ). Such variations are essentially due to modifications of the molecular 
structure of the PAHs, in different astrophysical environments. In particular, since each 
feature is attributed to a given vibrational mode, the ratio between these features will vary 
with quantities such as the charge, the hydrogenation or the size and shape of the molecule. 
The 3.3 jim PAH band arises from the radiative relaxation of CH stretching modes, while the 
11.3 /xm feature originates in the CH out-of-plane bending modes; CC stretching modes are 
responsible for the features between 6 and 9 /xm; CH in plane bending excitation produces 
part of the 8.6 /xm band. Now, laboratory studies and quantum calculations shows that the 
CC mod es are instrinsically weak in neutral PAHs, and become stronger when the PAHs ar e 
ionized dLanghoflj Il996l : kUamandola et al.l Il999l : iBauschlicherl I2OO2I : Ixim fc Savkallvl bood ) . 
Therefore, the 6 to 9 /xm bands will be much more intense for a PAH^ than for a PAH^, while 
it will be the opposite for the 3.3 and 11.3 /xm. Consequently, the ratios between the CC 
and the CH feature intensities depend on the charge of the PAHs, which is directly related 
to the physical conditions (e.g. intensity of the ionizing radiation field, electron density, etc.) 
in the environment where the emission is originating. 

Evidence of variations between features in different astroph ysical environments have 
been reported by many authors. For example, Ijoblin et al.l (119961 ) showed that the Is.e/Iii.a 
ratio decreases with increasing distance from the exciting star, in the refiect ion nebulae 



NGC 1333 - where Ia is the integrated intensity of the feature centered at A /xm. iHonv et al 



(I2OOII ) found a good correlation between the 3.3 and 11.3 /xm CH bands, in a sample of 
Galactic Hll regions, YSOs, and evolved stars, while they reported variations of I6.2/I11.3 

by a factor of 5^ The observations of Galactic and Magellanic Hll regions, presented by 

Vermeij et al.l (120021 ) . indicate that the ratios I6.2/I11.3, I7.7/I11.3 and Is.e/Iii.a are correlated. 
Furthermore, they suggest a segregation between the values of t hese ratios in the Milky 
Way and those in the Magellanic Clouds. iBregman fc Temil (120051 ) studied the variation of 
I7.7/I11.3 in three refiection nebulae. Assuming that this variation is controlled by the charge 
of the PAHs, they could relate this band ratio to the ratio Gy^lri p between the int e grate d 
intensity of the UV field, Go, and the electron density, rig. Similarly. ICompiegne et al.l (120071 ). 
studying the detailed variations of the mid-IR spectrum in the Horsehead nebula, attributed 
the high relative strength of the In. 3 feature to a high fraction of neu tral PAHs, due to the 
high ambient electron density. On the contrary, ISmith et al.l (120071 ) studied the variation 
of I7.7/I11.3 coming from the nuclear regions of the SINGS legacy program galaxies. They 
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find that this ratio is relatively constant among pure starbursts, but varies by a factor of 
5 among galaxies having a weak AGN. They interpret this effect as a selective destruction 
of the smallest PAHs by the hard radiation arising from the accretion disk, ruling out the 
explanation in terms of ionization of the molecules, in these particular environments. This 
interpretation is also supported by the high 11.3 /xm and the we akness of the 3.3 um band 
in the AKARI spectinm of the giant elliptical galaxy NGC 1316 (IKaneda et al.l 120071 ). 

The previous considerations stress the diversity of the possible interpretation of the mid- 
IR feature variations in galaxies. We need to identify the main physical processes controling 
the PAH bands, if we are to use them as diagnostic tools. This is the aim of this paper. It 
presents a quantitative analysis of mid-IR spectra {ISO and Spitzer) of Galactic regions, low- 
metallicity dwarf galaxies, quiescent spirals, starbursts and AGNs. To achieve our goal, we 
focus on identifying the main trends between the various PAH features, at different spatial 
scales and in different environments. We then link these variations to the ph ysical conditions 
insid e the studied region. Preliminary results of this study were published in iGalliand (12004 



20071). 



The paper is organized as follows. In ^ we present our sample and the data reduction. 
We discuss the spectral decomposition of the mid-IR spectrum, in ^ The results of this 
decomposition are presented in ^ we study the various trends between the band ratios 
within galaxies, and among different types of environment. Then, in ^ we provide a physical 
interpretation of these trends, when the structure of the ISM is resolved, and when it is not. 
Finally, we summarize our conclusions in §61 



OBSERVATIONS AND DATA REDUCTION 



2.1. The Sample 



In order to systematically study the properties of the mid-IR aromatic features, we 
include in our sample a wide variety of Galactic regions and galaxies, covering a large range of 
metallicities and star format i on act ivities. We merge the /a^O s ample s of starbursts and AGNs 
presented by Laurent et al.l (2000 ). spirals by Roussel et al. ( 2001 ). Magellanic regions and 
dwarfs bv iMadden et al.l kood ), and complement them by low-metallicity sources observed 
with Spitzer. We add to this sample, the ISO/ CAM spectra o f several Gal actic regions, like 
M 17 dCesarskv et aDll996bl). N GC 2023 dAbergel et al.lbood ). NGC 7027 dPersi et al.lllQQol ) 
and the Orion bar (ICesarskv et al.ll2000l). and the ISO/S WS spectra (2.5 — 45 /xm) of several 
compact Hll regi ons publishe d by IPeeters et al. (2002b). The SMC Bl#l spectrum is the 
one presented by iReach et al.l (120001 ). 
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The global properties of the selected sources are presented in Table [H If relevant, the 
distances were homogenised to Hq = 71 /i~^kms~^ Mpc~^. 
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Table 1. General properties of the sample. 
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R.A. 


Dec. 


Aperture 


Mid-IR 


Distance 


12 + 


log(0/H) 


Category 




(J2000) 


(J2000) 






spectrograph 


(Mpc) 




[ref.] 
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Table 1 — Continued 



Name 


R.A. 


Dec. 


Aperture 


Mid-IR 


Distance 


12 + 


log(0/H) 


Category 
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6" 


ISO/CAM 


4.5 


9.2 


m 


Spiral 


Circinus 


14^13^09^6 


-65°20'21'' 




40" 


ISO/CAM 


4.0 






SB/AGN 


Circinus e 


14^13^09^6 


-65°20'21'' 


40" 


-10" 


ISO/CAM 


4.0 






SB/AGN 


Arp 220 


15^34^57^2 


+23°30'11'' 




20" 


ISO/CAM 


73 






ULIRG 


Arp 220 p 


15^34^57^2 


+23°30'11'' 




10" 


ISO/CAM 


73 






ULIRG 


IR 15384 


15M2^17!1 


-53°58'3r' 


14" 


X 20" 


ISO/SWS 


2.7 kpc 






Hll region 


NGC 6240 


16^52^58!8 


+02°24'06'' 




20" 


ISO/CAM 


98 






SB/AGN 


NGC 6240 p 


16^52^58^8 


+02°24'06'' 




10" 


ISO/CAM 


98 






SB/AGN 


M 17 map 


18^20^22^0 


-16°12W 




6" 


ISO/CAM 


1.5 kpc 






PDR 


IR 18317 


18^34^24^9 


-07°54'4r' 


14" 


X 20" 


ISO/SWS 


4.9 kpc 






H II region 


NGC 6946 


20^34^51!2 


+60°09'ir' 




140" 


ISO/CAM 


5.5 


9.1 


[21] 


Spiral 


NGC 6946 p 


20^34^51^2 


+60°09'ir' 




40" 


ISO/CAM 


5.5 


9.1 


m 


Spiral 
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Table 1 — Continued 



Name 


R.A. 


Dec. 


Aperture 


Mid-IR 


Distance 


12 + log(0/H) 


Category 




(J2000) 


(J2000) 




spectrograph 


(Mpc) 


[ref.] 




NGC 6946 e 


20^34^51^2 


+60°09'ir' 


140'' - 40'' 


/50/CAM 


5.5 


9.1 [21] 


Spiral 


NGC 7027 


21^07^01^7 


+42°14'09'' 


14" X 20" 


ISO/SWS 


700 pc 




PDR 


IR 22308 


22^32^45^9 


+58°28'21'' 


14" X 20" 


ISO/SWS 


5.5 kpc 




Hll region 


IR 23030 


23^05^10^6 


+60°14'41'' 


14" X 20" 


ISO/SWS 


5.2 kpc 




Hll region 


IR 23133 


23^15^31^4 


+6i°oro8'' 


14" X 20" 


ISO/SWS 


5.5 kpc 




Hll region 


IR 23128 


23^15^46^0 


-59°03'ir' 


20" 


ISO/CAM 


180 




SB/AGN 


IR 23128 p 


23^15^46^0 


-59°03'ir' 


10" 


ISO/CAM 


180 




SB/AGN 


NGC 7714 


23^36^14^1 


+02°09'19'' 


10" 


Spitzer/IRS 


37 


8.5 [22] 


Dwarf 



References. 
[5] jAlloin et al 



- la 

(|l97i; 



Berg vail et 



P iHeckman et al 



al 



Otte et al 



Dufour et al 



Zaritskv et all (1994); 



(Il997h: [T3l iRavindranath &: PrabhiJ (|200lh: [ H 



200lh : [71 IPutil fc Rovl (Il999l): l8l Istorchi-Bergmann erall (1995h: 



([l998h: [TOl IRov WalshT (|l9^ ^ p ll IPilvugin et alJ (|20o4 : iKobulnickv Skill^ 



m 



m 



Mas-Hesse KunthI (Il999h : llTl llzotov Thuan 



Perez-Montero fc Diad (2003h: 



19981): Sch aerer et all hood): [T9l B resolin et all ([2OO4); 



Webster fc SmithI (|l983l ): [2T] iKobulnickv fc Johnsonl (|l999h : [22] IConzalez-Dekado et all (|l995[ ) 



Note. — The sources are ordered according to their right ascension. For ISO/ CAM sources, the apertures are 
circular, centered on the coordinates listed in the second and third columns. For ISO/SWS sources, the beam is 
14" X 20" (bands 1 and 2). The diameter of the aperture is given in the fourth column. When the diameter has the 
form Oi" — 62"^ it means that the we subtracted from the flux in the circular aperture of diameter ^1", the flux in 
a concentric circular aperture of diameter O2" < Oi" . The sources which are followed by "map" are those for which 
the signal-to-noise of the spectral-image is good; in this case, the aperture refers to the pixel field of view. The 
letter after the name of the source designates the type of aperture: when there is no letter, most of the emission of 
the object is encompassed in the aperture; "p" is for peak emission; "e" for extended emission. 
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2.2. ISO/CAm Data Reduction 



Most of the sources in Table [J were obs erved with ISO /CAM (ICesarskv et al.lll996al) on 



Madden et al. 



mm 



board the ISO sateUite (IKessler et al.lll996l ). We refer to the work of 
for a detailed description of the data reduction applied to all of the ISO/ CAM data used 
in this study. The CVF performed spectral imaging using a 32 x 32 detector array, with a 
sampling of 3''pixel~^ or 6^'pixel~^ in our cases, from A = 5 /xm to 16.5 /xm with one pointing 
of two CVFs, from A = 5 to 9.5 /xm and from A = 9.0 to 16.5 fim. The spectral resolution 
increases from A/ A A = 35 to 51 across the full spectra. 

For the ISO data treatment, we used the CAM Interactive Reduction (CIR, version 
AUGOl; IChaniall l2003l ). The subtraction of the dark currents was performed using the 
Biviano et al.l ( 19981 ) model which predicts the time evolution for each row of the detector, 
taking into account drifts along each orbit and each revo lution. We masked the glitches us- 
ing multi-resolution median filtering (IStarck et al.lll999l ) on each block of data after slicing 
the cube. Additional deglitching was performed manually, examining the temporal cut for 
each pixel. We corrected the systematic memory effects using the Fouks-Schubert method 
(jCoulais fc Abergell I2OOOI ) . We computed a hybrid flat-field image placing a mask on the 
source and computing a flat field outside this mask from the median of the temporal cut for 
each pixel. For the pixels which were on-source, the flat-field response was set to the cor- 
responding calibration flat-field. The conversion from Analog to Digital Units to mJy/pixel 
was performed using the standard in-flight calibration data base. To remove the sky contri- 
bution, sources smaller than the array were masked and, for a given wavelength, the median 
of the pixels which are off-source were subtracted from each pixel. For the more extended 
sources, we subtracted an independently observed zodiacal spectrum. The contribution of 
this spectrum was a free parameter varied in order to match the properly sky subtracted 
fluxes in the LW2 (6.7 fim) and LW3 (14.3 fim) broadbands. The final product is a 3D 
spectral-image. We integrated the spectrum using an aperture encompassing the entire 
galaxy, and obtained a ID spectrum. When the angular size of the source was larger than 
the array, we scaled the spectrum by matching IRASuj^m broadband flux with that derived 
from convolving the spectrum with the IRASi2,j^m bandpass. 



2.3. Spitzer/lUS Spectrum Extraction 



Several of the low-metallicity sources in Tabled] were not observed by ISO/CAM^ thus 
we complemented o ur database with mid-IR spectra fro ni the Syitzer/IRS spectrometer 
( Houck et al. 2004b ) on board the Spitzer Space Telescope ( Werner et al. 20041 ) . when these 
data were released. Among these galaxies, the spectrum of SBS 0335-052 has been published 
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bvlHouck et all fl2004al ). NGC 7714 by ISrandl et all (120041 ). NGC 5253 by iBeirao et al 
( 120061 ). Haro 11 and I Zw 18 bv IWu et al.l (120061 ). We refer to these studies for a detailed 
presentation of their mid-IR properties. We considered only low-resolution data, taken with 
the SL (Short-Low) module, from A = 5.2 /xm to 14.5 /xm, and the LL (Long-Low) module, 
from A = 14.0 fim to 38.0 /xm, both with a spectral resolution of A/AA 64 — 128. 

We retrieved the Basic Calibrated Data (BCD) that have been preprocessed by the 
Spitzer Science Center (SSC) data reduction pipeline. They have been converted to flux 
density, corrected for stray light and flatflelded. The extraction of the spectra from the 2D 
space/ wavelength images was performed with the Spectral Modeling, Analysis and Reduction 



Tool (SMART, version 5.5.6: iHigdon et al.ll2004l ). We first inspected the BCD images and 
identified the hot pixels which had not been masked out by the SSC. We replaced them by the 
median of their neighbors. For each module, one order is on-source and the other is off-source, 
at a time. Then the positions are switched. We subtracted the off-source spectrum from 
that of the source, in order to remove the sky emission. The extraction of the ID spectrum 
was performed inside a column whose width varies with wavelength. We have excluded the 
bonus order. Then, the various frames, for each nod position, were coadded. Since the long 
wavelength end of the SL module and the short wavelength end of the LL module were not 
systematically overlapping, we finally scaled the SL module, in order to obtain a continuous 
spectrum. This scaling factor can be as large as 50%. At the time of this publication, 
the Spitzer/IRS data handbook recommends that the signal to noise ratio not be derived 
from the uncertainties generated by the pipeline. Instead, we adopted the recommended 
systematic error of 20%. To take into account these statistical variations, we smoothed the 
spectra into a AA 0.2 /xm window (4 points), and took the standard deviation inside this 
window as the error. As we did with the ISO/ CAM spectra ( §2.21) . in order to compensate 
for the fact that we may be overlooking some extended emission, we scaled the spectrum to 
match the IRASu/^m and IRAS^.^^.r n . In the case o f II Zw 40, the Spitzer /lYlS spectrum does 
not exhibit the PAH features that iMadden et al.l (120061 ) detected in the extended emission, 
since the Spitzer/lYl^ observations were conducted in staring mode. Indeed, this Spzteer/IRS 
observation does not encompass the region where the tail was detected. Thus we prefer to 
use the ISO/ CAM data for this galaxy. 



3. THE SPECTRAL DECOMPOSITION METHODS 



3.1. Inventory of the Physical Components 



We restrain our study to the spectral range 5 — 16 /xm, since it is the overlap between 
the various data sets we have compiled. The main physical components contributing at these 
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wavelengths are the foUowing. 

• The ionic gas emits strong fine-structure ionic fines. Tfie [Nellljis.se^m, [Nell]i2.8i/im5 
[Siv]io.5i^m, [Arlll]8.99Mm and [Arll]6.98Mm are tfie brigfitest fines. 

• Several ubiquitous broad features, originating from tfie molecular modes of stocfiasti- 
cally fieated PAH molecules. Tfie most prominent of tfiem are centered at 6.2, 7.7, 8.6, 
11.3 and 12.7 /xm. Tfie [Ne Il]i2.8i/im line is blended witfi tfie 12.7 /xm feature. 

• Very small grains (VSG), fluctuating in temperature around a few fiundred degrees 
Kelvin, produce a continuum emission. In addition, fiot grains in H ll regions, at 
tfiermal equilibrium witfi tfie radiation field, may contribute to tfiis continuum. 

• Tfie extinction feature at 9.7 /xm, attributed to amorpfious silicates, can be seen toward 
tfie most embedded sources. 

• Molecular lines are also present, especially tfie H2 0-0 S(3)9.7^m line. However, tfiey are 
mucfi weaker tfian tfie otfier components. 

Fig.[T]sfiows tfie variations of tfie aromatic feature spectrum among galaxies (top panel), 
inside one galaxy (middle panel) and from a tfieoretical point of view (bottom panel). Tfiis 
figure sfiows tfiat tfie 6.2, 7.7 and 8.6 fim features are qualitatively tied togetfier, and tfiat 
most of tfie variations in tfie mid-IR aromatic spectrum is a variation of tfiese tfiree features 
relative to tfie 11.3 /xm band. Tfie lower panel of Fig. [1] suggests tfiat most of tfiese variations 
could be explained by a variation of tfie neutral-to-cationic-PAH ratio. In wfiat follows, we 
quantify tfiese variations in our sample, in order to explore tfie validity of tfiis explanation. 



3.2. Measuring the Intensity of the Aromatic Features 



Measuring tfie intensity of tfie aromatic features is uncertain, due to tfie intrinsic widtfi 
(AA 0^ 1 /xm) and tfie complexity of tfie band profiles. Indeed, contrary to gas lines, a 
large fraction of tfie energy of tfie PAH bands is radiated in tfie wings. Tfiese wings extend 
far outside of tfie central wavelengtfi and can be difficult to reliably disentangle from tfie 
underlying continuum emission. Tfierefore, a proper extraction of tfie feature requires an 
assumed band profile. However, tfie actual profi le of eacfi band is not known. It is asymmetric 



and varies as a function of environment (e.g. IPeeters et al.l l2002al : Ivan Diedenfioven et al 



2004 ) . Consequently, several profile s fiave been used i n tfie literature, motivated by different 
pfiysical arguments. For example, iBoulanger et al.l (119981 ) proposed tfiat tfie widtfi of tfie 
band is a consequence of tfie continuous redistribution of intramolecular vibrational energy. 
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[Aril] [Arm] [SIV] [Nell] [Nelll] 




4 6 8 10 12 14 16 
Rest frame wavelength [/xm] 



Fig. 1. — Various spectra normalized to the 11.3 /xm feature. The top panel shows the 
total integrated spectra of a low-metallicity galaxy (UM 448), a normal spiral (M 51), and a 
LIRG (NGC 4945). The middle panel shows the spectra of three different regions inside the 
starburst galaxy M 82. The bottom panel shows, for comp arison, the absorption coe fficient 
of neutral and cationic PAHs measured in laboratory by lAUamandola et al.l ( 119991 ) . The 
solid vertical lines mark the wavelengths of the brightest lines, and the dashed vertical lines 
mark the wavelengths of the major bands. In this figure and in what follows, Fjy{X) is the 
monochromatic flux density at wavelength A /xm. 
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between different excitation levels, at high temperature. They estimated the transition 
timescales to be < 10~^^ s, and showed that the PAH ban ds were well repre sented by 
Lorentzian profiles. Conversely, a Drude profile was used bv iLi fc Draind (120011 ). Such a 
profile describes the electric permeability in a solid or a large molecule. Finally, numerous 
studies simply measure the tip of the aromatic band, and assumes that the total fiux radiated 
in the band scales with this quantity. 



Spline method 

[Aril] [Arm] [SIV] [Nell] [Nelll] 




400 



300 



200 



4 6 8 10 12 14 16 
Rest frame wavelength [/xm] 



Lorentzian method 

[Aril] [Arm] [SIV] [Nell] [Nelll] 




400 



300 



200 



4 6 8 10 12 14 li 
Rest frame wavelength [/xm] 



Fig. 2. — Demonstration of the two spectral decomposition methods for the total spectrum 
of M 82 and 30 Doradus. 



In order to take into account the uncertainty of the actual shape of the aromatic bands, 
we systematically analysed our spectra using two different spectral decomposition methods. 
These two methods differ in the way the continuum and the aromatic bands are fitted. 
They are demonstrated in Fig. [21 on a PAH-rich spectrum (M 82), and PAH-poor one 
(30 Doradus). 



The Syli ne method (hereafter identified by the letter 5), used e.g. by iHony et al.l ( 120011 ) 
and IVermeij et al.l ( 120021 ) , takes into account only the tip of the aromatic bands when 
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computing their intensities, as shown in left panels of Fig. [21 The continuum (grey 
filled curve) is fitted with a spline function, constrained at the wavelengths A = 5.04, 
5.47, 5.84, 9.18, 10.85, 13.82, 14.60, 15.08, 16.00, 16.15 //m. A second spline function 
is fitted to the spectrum, corresponding to the previous continuum, plus the "p/aieaix" 
under the 6.2, 7.7 and 8.6 fim bands, and the plateau under the 11.3 and 12.7 /xm 
bands. In addition to the previous wavelengths, this spline function is also constrained 
at A = 6.61, 6.80, 7.13, 8.24, 8.80, 11.85, 12.20, 13.22 //m. The diflFerence between the 
total spectrum and this continuum+plateau component defines the lines and features 
(line filled areas of the left panels of Fig. [2]). This method is relatively robust. However, 
it does not allow us to separate the [Ne Il]i2.8i/im from the 12.7 fim band (most of our 
spectra have a low spectral resolution). 



The Lorent z ian method (hereafte r iden tified by the lett e r >C), used e.g. bv iBoulanger et al 
(119981 ) , iLaurent et al.l (120001 ) and IVerstraete et al.l (120011 ) , assumes that the aromatic 
bands have lorentz profiles. The wings of these lorentz ians account for th e underlying 
plateau. A variation of this method, implemented by ISmith et al.l (120071 ). consists of 
modeling the PAH features with Drude profiles. We proceed as following. 

1. The [Nelll]i5.56^m, [Nell]i2.8i/im, [Siv]io.5i^m, [Arlll]8.99Mm, and [Arll]6.98/im ionic 
lines are fitted with gaussian functions (diagonal line filling in the right panels of 
Fig.[2j). For each (number i) of the N^^^ lines, the central frequency, lyf^^^ is fixed, 
the width is fixed by the spectral resolution of the instrument, Az/gpectro, and the 
total fiux of each line, F-^^^^ is the free parameter. 

2. The PAH bands are fitted with lorentzian functions (horizontal line filling in the 
right panels of Fig. [2j). The considered features are centered at: A = 5.3, 5.7, 6.2, 
7.7, 8.6, 11.3, 12.0, 12.7, 13.6 and 14.3 //m. For each (number i) of the A^band 
band, the central frequency, i/^^^^^ the width, Az/]^^^"^, and the total fiux, F^^^^^ 
are the free parameters. We emphasize the fact that we fit only one band around 
7.7 /im, where there are actually two bands centered at 7.6 /xm and 7.8 /xm. We 
proceed this way, in order to keep the number of parameters reasonably low. We 
will use the variation of the centroid of the 7.7 fim feature, noted A7.7, as an 
indicator of the 7.6 fim to 7.8 fim ratio. 

3. The VSG continuum is fitted with the sum of Acont modified black-bodies having 
the absorption efficiency, Qahsi^)^ of gra phite in the RayleiRh approximation (grey 
filled area on the right panels of Fig. [2{ iLaor fc Draindll993l ). The temperatures, 
T/^^^*, and the total fiuxes, F^^^^^ are the free parameters. We adopt Acont = 2. 

The suni of th e previous components is multiplied by exp(— r), where rjy ) is the 



MathisI (jl990l ) extinction law, with the iDudley fc Wynn-Williamsl ( 119971 ) silicate 
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features. The column density is the only free parameter. 
All these components are fitted simultaneously as part of one mathematical function: 
F,{u) = {Fl'-'iiy) + F^^-'^iu) + Fr\v)) x exp {-t{v)) , (1) 

with: 



i=l 



line 



^^^spectro 



exp —2 



[jy — jyl 



line A 2 



A 7/2 

spectro 



^band -i 

^band 



,band 



i=l 



A^cont 



;^cont 



cont 



(2) 



J c/ (16 iJ,m) 

the various Fjy being the monochromatic flux densities. An advantage of the Lorentzian 
method over the Spline method is that it aUows us to separate the [Nell]i2.8i/im line and 
the 12.7 fim feature, and to study the variations of the centroids of the bands. However, 
when the PAH-to-VSG ratio is very low, as in the case of 30 Doradus (Fig. [2j), the 
width of the features is more uncertain. 



In both cases, Icont is defined as the integrated intensity of the continuum between 10 and 
16 /im, and Ipah, the sum of the intensities of aU the bands between 5 and 16 /xm. We 
emphasize that our methods work automaticaUy, without any by eye adjustment. This 
condition prevents systematic effects that could result from arbitrary ch oices of parameters . 



These two methods are orthogonal but give similar trends (see ^ and ISmith et al.l 120071 ). 
The Spline method is clearly continuum-biased in its assumptions whereas the Lorentzian 
method is a line-biased assumption method. 

Our two methods can be applied to an integrated spectrum, or to each pixel of a spectro- 
image, provided that the signal-to-noise ratio is suflicient. We have systematically applied the 
two methods to all the spectra presented in Table [H In our analysis, we will systematically 
consider the results of the two methods before drawing conclusions on the PAH properties. 
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4. SYSTEMATIC ANALYSIS OF THE MAJOR BAND RATIOS 

4.1. Correlations Exhibited Among Integrated Spectra 

Figs. [3H11 show select correlations between band ratios of the integrated spectra of our 
sample (Table [1]), obtained with the two methods presented in §3.2[ We focus on the four 
brightest bands at 6.2, 7.7, 8.6 and 11.3 /xm, and study the correlations between the various 
ratios. The intensities of the features are reported in App. [A] (Table Hj), and the parameters 
of the various correlations are given in Table [21 As mentioned in §3.21 the two methods are 
uncertain for very low values of the PAH-to-VSG ratio. Thus we define two subgroups of 
data. 

1. The data that we consider to be reliable are marked with black symbols in Figs. [HUH 
We define them as the measurements which have a signal-to-noise ratio, at A = 7.7 /xm, 
larger than 6, and IpAn/Icont ^ 0.5 with the Lorentzian method^ and IpAn/Icont ^ 0.3 
with the Spline method. 

2. The fits that we consider to be less certain are the complementary data (grey symbols 
in Figs.[3Hll), having lower signal-to-noise and PAH-to-VSG ratios. 

Fig. [3] shows that there is an excellent linear correlation between the ratios I6.2/I11.3, 
I7.7/I11.3 and Is.e/Iii.s, using either of the two methods. The variations of these ratios spread 
roughly over one order of magnitude, while the ratios I6.2/I7.7, and I7.7/I8.6 are roughly 
constant within the error bars (Fig. H]). In general, the measure of the intensity of the 

8.6 fim band is less accurate than for the other main bands. Indeed, it is less intense than 
the 6.2 and 7.7 /xm features. In addition, it is merged with the long wavelength wing of the 

7.7 /im band. Finally, this particular feature is significantly affected by the silicate extinction 
feature around 9.7 /xm, when the source is deeply embedded. This extinction is not corrected 
in the case of the Spline method^ and is corrected very simply in the case of the Lorentzian 
method. That is the reason why correlations involving the 8.6 /xm feature are always more 
dispersed than the others. Quantitatively, the correlation coefficients are around 0.8 for the 
four panels of Fig. [3l except for that involving the Is.e feature, with the Spline method. Most 
of the outsiders are based on fits that we consider to be uncertain. 

Fig. [5] shows the consistency between the two methods. It demonstrates that, whatever 
method we use to measure the band ratios, the order of the various measures is conserved; 
i.e. a spectrum A having a lower I6.2/I11.3 ratio than a spectrum i3, using the Spline method^ 
will also have a lower ratio than B using the Lorentzian method. Moreover, the fact that 
both methods give similar trends means that we did not artificially incorporate part of the 
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continuum flux within the aromatic band intensities. Indeed, it is possible to accidentaUy 
account for a fraction of the continuum intensity within the wings of the PAH profile, with 
the Lorentzian method^ especially within the 11.3 /xm feature. However, this bias is not 
possible with the Spline method^ since it integrates only the tip of the band. Therefore, the 
agreement between the two methods allow us to rule out this bias and claim that the trends 
of Fig. [3] are not induced by the fitting methods. 

These first relations, established on integrated spectra, indicate that the properties 
of the PAHs throughout different types of galaxies and Galactic regions are remarkably 
homogeneous. They are consistent with a significant variation of the 6.2, 7.7 and 8.6 /xm 
features relative to the 11.3 /xm band, coupled with an absence of significant variations 
among the 6.2, 7.7 and 8.6 /xm bands. 

Finally, as a consistency check. Fig. [6] shows the variations of the band ratios as a 
function of distance. The correlation coeflicient is 0.23 for the Lorentzian method^ and 0.20 
for the Spline method. This absence of correlation shows the uniformity of the sample and 
rules out the possibility of variations induced by aperture eflFects. 

4.2. Spatial Variations of the Band Properties 

We now study the spatial variations of the properties of the mid-IR features, by analyzing 
the spectral maps of a sub-sample of sources: IC 342, M 17, M 51, M 82, M 83, 30 Doradus, 
and the Orion bar. These objects are those which are spatially resolved, satisfy the two 
PAH-to-continuum and signal-to-noise ratio thresholds listed in §4.11 and show significant 
variations of the band ratios. Similar to §4.11 we systematically compare the results of the 
two methods, but instead of applying them to integrated spectra, we fit the spectrum of 
each pixel of the spectral maps. We then degrade the images of each spectral component to 
the spatial resolution of the longest wavelength of the spectrum (FWHM 9^'). We achieve 
this by convolving the image of a given component (e.g. I6.2), with the PSF at the longest 
wavelength (i.e. A = 16 /xm) deconvolved beforehand by the PSF of the component (i.e. 
A = 6.2 /xm). Figs. El to [12] (as weU as Figs. [22] to [36] in App. \E) show the same correlations 
as in Figs. [HHSl obtained inside resolved sources, instead of integrated spectra. Each data 
point represents the fit of the spectrum of one pixel. We selected the points, according to 
the criterion defined in §4. It the signal-to-noise ratio at A = 7.7 /xm must be larger than 
6, and IpAn/Icont ^ 0.5 with the Lorentzian method^ and IpAn/Icont ^ 0.3 with the Spline 
method. The pixels which do not fall into this category have not been considered here. 
The parameters and statistical quantities relative to these correlations are given in Table [21 
Figs. [13] to [H] (as well as Figs. [37] to [H] in App. \^ show the spatial distributions of the 



-19- 



Lorentzian Method 



O Dwarfs 
^ Spiral 
^ SB/AGN 
□ HII/PDR 





0.1 



Spline Method 




0.1 



Fig. 3. — PAH properties for integrated spectra. The black error bars are for the mea- 
surements that we consider to be reliable, while the grey are more uncertain. In each 
panel representing a relation Y versus X, the grey filled area is the domain between 
Y = ayjx X X — IcTy/x, and Y = ay/x x X + Idy/Xi Y = ay/x X X being the linear 
correlation of the data, and ay/x^ the dispersion of the data around this correlation. 
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Lorentzian method Spline method 




Fig. 4. — PAH properties for integrated spectra (continued). The same symbol conventions 
are adopted as in Fig. [31 
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^ea/^iis (Lorentzian) 

Fig. 5. — Comparison between the two methods. The same symbol conventions are adopted 
as in Fig. [3l 
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Fig. 6. — Correlation of the band ratio with the distance. The absence of correlation shows 
the uniformity of the sample. The same symbol conventions are adopted as in Fig. [31 

components. 

In general, the correlations presented in Figs. [3H5] still hold inside individual objects 
as well. Therefore, the variations of the PAH band ratios are independent of the spatial 
resolution. For example, the pixel size is c:::^ 0.1 pc in M 17 and the Orion bar, while it is 
1 kpc in external galaxies. From the images, we notice that the ratio I6.2/I11.3 (as well as 
the ratios I7.7/I11.3 and Is.e/Iii.a) is roughly correlated with the intensity of the PAH emission, 
the highest ratios being found in the brightest regions (spiral arms, starburst region, etc.), 
while the lowest ratios are generally found in the extended, low-luminosity regions. We do 
not detect any significant variations of the centroid of the 7.7 /xm feature inside each source. 
However, the average centroid varies from one source to another. 

App. El comments on each individual source. Here we will illustrate the systematic 
variations within the individual sources by focusing on two extragalactic sources: the edge- 
on irregular starburst galaxy M 82 and the face-on spiral galaxy M 51. For M 82 (Fig. [T3l) . 
the maximum values of the I7.7/I11.3 ratio are found along the disc, in the star forming region, 
where the infrared emission is the highest. The value of the band ratio drops by a factor of 
3 at 200 pc above and below the disc, in the halo. For M 51 (Fig. [HI), the I7.7/I11.3 ratio 
is maximum along the circumnuclear ring. It is somewhat lower in the 10'^ nuclear region. 
The band ratio is also high along the spiral arms, coinciding with the infrared bright blobs, 
but can exhibit variations by a factor of 2. Finally, the ratio drops in the interarm region. 
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The highest values of the I7.7/I11.3 — 7 ratio {Lorentzian method)^ found in the nuclear 
star forming regions of M 82 and M 51, are similar to the values found in the compact 
Galactic Hll regions IR 15384 and IR 18317 (Tabled]). Intermediate values of I7.7/I11.3 — 3 
(Lorentzian method)^ similar to the Galactic reflection nebula NGC 2023 are found in the 
halo of M 82 and along spiral arms of M 51. The lowest values of I7.7/I11.3 ^ 1 {Lorentzian 
method)^ seen in the planetary nebula NGC 7027, correspond to the outermost regions of 
the two galaxies. 
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M 82 (Lorentzian method) 




M 82 (Spline method) 




Fig. 7. — PAH band ratios within M 82. The grey filled areas are the correlations obtained 
for the integrated spectra (Fig. [3]). The correlations obtained inside M 82 are represented 
by the hatched region. The white symbol is the value of the global measurement over the 
entire galaxy. 
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Fig. 9. — Comparison between the two methods in M 82. The same symbol conventions are 
adopted as in Fig. [71 
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M 51 (Lorentzian method) 




Fig. 10. — PAH band ratios within M 51. The same symbol conventions are adopted as in 

Fig.m 
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Fig. 12. — Comparison between the two methods in M 51. The same symbol conventions 
are adopted as in Fig. [71 
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M 82: Ip.Htot & Ip.h(6.3)/Uh(11-3) spime M 83: & Ip,H(6.a)/Ip,H( 1 1 -3) 

I \ \ \ \ 1 I \ \ \ 




Fig. 13. — Spatial distribution of the PAHs in M 82. For each method, the image is the 
total PAH intensity, and the contours are the I6.2/I11.3 ratio. The shaded circle indicates the 
beam size. 
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Fig. 14. — Spatial distribution of the PAHs in M 51. For each method, the image is the 
total PAH intensity, and the contours are the I6.2/I11.3 ratio. The shaded circle indicates the 
beam size. 



Table 2. Parameters of the various correlations. 



Source 


X 


Y 


Method 






(^Y/X 


PY/X 


Integrated spectra 


I6.2/I11.3 


I7.7/I11.3 


S 


1.78 





0.69 


0.89 


(Fig.© 






C 


2.70 





0.64 


0.85 




Is.e/Iii.s 


I7.7/I11.3 


S 


8.21 





1.09 


0.44 








C 


4.43 





0.79 


0.78 




I6.2/I11.3 


I6.2/I11.3 


s/c 


1.03 





0.46 


0.45 


M 82 


I6.2/I11.3 


I7.7/I11.3 


s 


1.86 





0.32 


0.97 


(Fig. [71) 






c 


2.70 





0.32 


0.97 




I8.6/I11.3 


I7.7/I11.3 


s 


8.06 





0.53 


0.87 








c 


4.04 





0.39 


0.95 




I6.2/I11.3 


I6.2/I11.3 


s/c 


1.07 





0.23 


0.88 




I7.7/I11.3 


A7.7 [iim] 


c 





7.68 


0.012 




IC 342 


I6.2/I11.3 


I7.7/I11.3 


s 


1.50 





0.63 


0.82 


(Fig. [22]) 






c 


2.33 





1.35 


0.79 




I8.6/I11.3 


I7.7/I11.3 


s 


5.71 





0.69 


0.81 








c 


1.89 





1.41 


0.66 




I6.2/I11.3 


I6.2/I11.3 


s/c 


1.12 





0.40 


0.74 




I7.7/I11.3 


A7.7 [l-im] 


c 





7.68 


0.012 




M 51 


I6.2/I11.3 


I7.7/I11.3 


s 


1.49 





0.50 


0.62 


(Fig.nni) 






c 


2.59 





0.58 


0.65 




I8.6/I11.3 


I7.7/I11.3 


s 


7.19 





0.62 


0.41 








c 


4.55 





0.56 


0.65 




I6.2/I11.3 


I6.2/I11.3 


s/c 


1.01 





0.42 


0.33 




I7.7/I11.3 


A7.7 [l^rn] 


c 





7.72 


0.012 




M83 


I6.2/I11.3 


I7.7/I11.3 


s 


1.30 





0.52 


0.62 


(Fig. ESI 






c 


2.19 





0.54 


0.67 




I8.6/I11.3 


I7.7/I11.3 


s 


7.18 





0.69 


0.43 








c 


3.52 





0.70 


0.66 




I6.2/I11.3 


I6.2/I11.3 


s/c 


1.07 





0.41 


0.39 



Table 2 — Continued 



Source 


X 


Y 


Method 






(^Y/X 


PY/X 




I7.7/I11.3 


A7.7 [l-im] 


C 





7.71 


0.016 




30 Doradus 


I6.2/I11.3 


I7.7/I11.3 


S 


1.43 





0.79 


0.78 


(Fig.EHl) 






C 


2.31 





0.67 


0.68 




Is.e/Iii.a 


I7.7/I11.3 


S 


11.48 





1.61 


0.35 








C 


2.10 





0.62 


0.49 




I6.2/I11.3 


I6.2/I11.3 


s/c 


2.01 





0.86 


0.44 




I7.7/I11.3 


A7.7 [iim] 


c 





7.66 


0.016 




M 17 


I6.2/I11.3 


I7.7/I11.3 


s 


1.96 





0.33 


0.81 


(Fig. El 






c 


2.41 





0.45 


0.93 




I8.6/I11.3 


I7.7/I11.3 


s 


11.91 





0.81 


0.33 








L 


4.U4 





r\ A 

0.84 


0.80 




I6.2/I11.3 


I6.2/I11.3 


S/C 


1.01 





0.25 


0.61 




I7.7/I11.3 


A7.7 [m^] 


c 





7.69 


0.012 




Orion bar 


I6.2/I11.3 


I7.7/I11.3 


s 


1.77 





0.32 


0.83 


(Fig.[3lD 






c 


2.32 





0.21 


0.95 




I8.6/I11.3 


I7.7/I11.3 


s 


8.38 





0.38 


0.91 








c 


3.72 





0.52 


0.86 




I6.2/I11.3 


I6.2/I11.3 


s/c 


1.01 





0.29 


0.40 




I7.7/I11.3 


A7.7 [iim] 


c 





7.68 


0.011 





Note. — This table gives, for two given measured quantities X and F, the 
parameters corresponding to the fit oiY hy ay/x x X + hy/x- The dispersion of 
the measures around this correlation is (Jy/x = \/ ((^ — (^Yjx x X + ^y/x)^)- 
In the case where hyjx =0, pyjx is the Unear correlation coefficient. 
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INTERPRETATION OF THE AROMATIC FEATURE VARIATIONS 



5.1. Origin of the Band Ratio Variations 

The main conclusions of §1] are that (i) the ratios I6.2/I11.3, I7.7/I11.3 and Is.e/Iii.a are 
correlated and span one order of magnitude throughout our sample, and (ii) the ratios 
I7.7/I6.2 and I8.6/I6.2 are not significantly anticorrelated with the I7.7/I11.3 ratio, and do not 
show significant variations within the sources which have the best signal-to-noise ratios 
(especially M 82). As explained in and demonstrated in Fig. [U the intensity of the 
features between 6 and 9 /xm relative to the 11.3 jim band are one order of magnitude higher 
for PAH+ than for PAH°. Therefore, the universal correlations displayed from Fig. [3] to 
Fig. [33] can be attributed to variations of the charge of the carriers of the aromatic features. 
However, before concluding, we first need to explore the ability of other physical processes 
to reproduce these trends. 

First, dehydrogenation of the PAHs have an eff ect similar to ionization, on the mid- 



IR spectrum, as shown by several laboratory studie s (lAUamanc 



19961 : Ijochims et al.lll999l ). However, as reported bv lAUain et al. 



ola et 



al 



19891 : 



AUain et al, 



fll996h . only PAHs contain- 



ing less than Nq — 50 carbon atoms can be considerably dehydrogenated. At the same time, 
this threshold corr esponds to the minimum size of PAHs that can survive in most PDRs 



( AUain et al.l 1 19961 ). For comparison, the PAH size distribution derived by IZubko et al. 

for the diffuse Galactic ISM, has a lower cut-off of Nq^^ = 20. Their emission is 
dominated by PAHs of Ac 50 — 500, depending on their excitation rate. This conclusion 
is supported by the observed linear c orrelation of the 3.3 fim band with the 11.3 //m band 
in Galactic sources (IHony et al.ll200ll ). If dehydrogenation were important then a non-linear 
relation would be expected as duo and trio gro ups (with bands lonRwards of 11.3 um ) were 
converted into solo groups (with a 11.3 /xm band: ISchutte et al .111 9931 : iHonv et al. Like- 
wise, if dehydrogenation were important, then the decrease in the I11.3/I6.2 ratio (decreasing 
H coverage) would be expected to be accompanied by a decr easing Ii^ r/In .s r atio (conver- 
sion of duo/trio's into solo's) in contrast to the observations (IHony et al.ll200ll ). Therefore, 
we can neglect the effect of dehydrogenation on our mid-IR spectra. 

Second, deep extin ction by the sil i cate fea t ure at 9.7 /ini can ca use variations of t h e ban d 
ratios, as proposed by IPeeters et al.l (j2002bl ). ISpoon et al.l (120021 ) or iBrandl et al.l (120061 ). 
Indeed, the 11.3 /xm feature is more absorbed than the 7.7 fim band. This effect would 
also cause the 8.6 /xm feature to suffer from the same amount of extinction as the 11.3 /xm 
band, and these two bands should be correlated. However, this is not what Figs. [3] to Fig. [35] 
show. Moreover, the Lorentzian method corrects the bands for extinction (§3j). The good 
agreement between the Lorentzian method and the Spline method tells us that extinction can 
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not explain the majority of the observed band ratio variations. FinaUy, there is no evidence 
for sihcate absorption in the spectra of the Orion bar. Indeed, the total gas and dust column 
along the line of sight for this source is insufficient to produce noticeable 10 /xm absorption. 



3 
2.5 



Galactic 



Young burst 



Nc =10; XisRF=l 




8 10 12 
Wavelength [/xm] 



14 



8 10 12 14 

Wavelength [/xm] 



Fig. 15. — Effects of the temperature distribution on the PAH spectra. The various curves 
are the theoretical emission spectra, varying the shape (Galactic or young starburst) and 
the intensity of the ISRF (xisrf), as well as the lower cut-off of the size distribution (N^^^). 
The size distri bution is the BARE -GR-S of IZubko et al.l (120041 ). The absorption efficiencies 
are taken from lDraine fc Lil (120071 ). 



Third, a modification of the temperature distribution of the PAHs can affect the relative 
weight of the aromatic features. For example, truncating the highest temperatures would 
have the effect of decreasing the intensity of the short wavelength bands relative to the long 
wavelength ones. There are essentially two ways of modifying the temperature distribution 
of the molecules: (i) hy varying the interstellar radiation field (ISRF) intensity or shape, 
and (ii) by changing their size distribution. The latter effect could happen in the vicinity 
of a very hard radiation source, like an AGN, whe re the smallest PAHs could undergo 
photosublimation, as proposed by ISmith et al.l (120071 ). 



We performed theoretical modeling of the PAH emission spectra, in order to study 
the latter effects. These spectra, which take into account the stocha stic heating of the 
particles, are shown in Fig. [151 We adopted the PAH size distribution bv lZubko et al.l (120041 . 
bare grains, solar abundance constraints; BARE-GR-S), and the absorption efficiencies by 
Draine fc Lil (120071 ). with an ionization fraction of 50%. We consider successively two classes 
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of ISRF, in order to expl ore the effect of the h ardness of the radiation on the PAH emission: 
(i) the Galactic ISRF of iMathis et al.l (119831 ) and (ii) d. very hard ISRF, corresponding to 
an instantaneous burst of st ar formation with a Salpeter IMF , synthesized with the stellar 
evolutionary model PEGASE (IFioc fc Rocca-Volmerangd 119971 ). The latter is normalised to 
the intensity of the iMathis et al.l (119831 ) ISRF. We also vary the intensity of each ISRF, 
by multiplying them by a factor Xisrf = 1 — 10^. In addition, we explore the effect of a 
possible small PAH segregation, by varying the value of the lower cut -off of the PAH size 



distribution. This lower cut-off is N^^^ = 20 carbon atoms, for the IZubko et al.l (120041 ) 
model. We vary it up to A^^?^^ = 10^. Finally, we apply a screen extinc tion to the n i odeled 
spectra by multiplying them by a factor exp[— r(A, ^y)]. We adopt the IZubko et al.l (120041 ) 
optical depth, r(A, Ay), and plot the results using two values of the V band attenuation: 
Ay = (no extinction) and Ay = 10. 



1 1 1 

- N^^''=20; Mathis et al. 


(1983) 


■ __. j^min^^Q3. y^g^^YllS et al. 


(1983) [ 


- N™"'=20; Young burst 




- - - ■ Nf''= 10^; Young burst 














1 1 1 


1 



1 10 100 1000 10^ 10^ 



XlSRF 



Fig. 16. — Theoretical modeling of the effect of the PAH temperature distribution on their 
band ratios, with a fixed ionization fraction (50%). No extinction is applied. 

Figs. fTGlfTTl synthesize the results of this modeling. Fig. [T6l demonstrates that effect of the 
radiation field intensity on the I7.7/I11.3 band ratio, for the two types of radiation fields and 
the two extreme PAH size cut-offs. This figure shows that, for a given size distribution and 
ISRF type, the I7.7/I11.3 ratio is independent of the intensity of the ISRF, up to Xisrf — 10^. 
Indeed, the PAHs are stochastically heated below Xisrf — 10^, therefore the shape of their 
emission spectrum is independent of the ISRF. However, when Xisrf ^ 10^, the largest 
PAHs reach thermal equilibrium, and the temperature fluctuation range shrinks significantly 
enough to affect the ratio between the bands. Fig. [TTl shows how every effect considered here 
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10 




Fig. 17. — The two panels present the results of the same modeling as in Fig. [161 the only 
difference being that the top panel shows the resulting I7.7/I6.2, while the bottom panel shows 
the I8.6/I6.2- To read this figure, first look at the black empty polygon with the circled corners. 
This polygon shows the effect of varying both the minimum cut-off PAH size, A^^?^^, and the 
ISRF intensity, Xisrf, with the Galactic ISRF (labeled Galaxy) and with no extinction. The 
values between parenthesis, close to each corners are the values of these two parameters at 
the corner. Second, the grey empty polygon is the analog of the previous one, except that 
the ISRF is now the Starburst one (labeled Burst). Finally, the two grey filled polygons 
are the analog of the two previous ones, except that we applied a screen extinction with 
Ay = 10. As a summary, the various arrows show the sense of variation of the band ratios 
with each parameter, including the ionization degree of the PAHs. 
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affects the relation between I7.7/I11.3 and I7.7/I6.2 and I8.6/I6.2- The arrows indicate the 
direction toward which these parameters affect the relation. This figure shows that the 
combination of all the effects affecting the PAH temperature distribution (ISRF intensity, 
ISRF hardness and size distribution) can not explain a variation of the I7.7/I11.3 ratio by 
more than a factor 3. Moreover, it shows that if these effects were responsible for the 
variation of the I7.7/I11.3 ratio, then the I7.7/I6.2 ratio would also be anticorrelated with it. 
The ratio I7.7/I6.2 is particularly very sensitive to A^™. According to Fig. [TTl an increase of 
the I7.7/I11.3 ratio by a factor 3 would imply a decrease of the I7.7/I6.2 ratio by a factor of 
2. 

One of the strong results of §4.11 and §4.21 is that the ratios I7.7/I6.2 do not show any 
trend with I7.7/I11.3. Therefore, the correlations discussed in §l]can not be attributed to a 
modification of the temperature distribution. In particular, we can rule out an interpretation 
of the variation of the I7.7/I11.3 in term s of destruct i on of the smallest PAHs. These results 
are in good agreement with those by iHony et al.l ( 120011 ). who studied a wide sample of 
Galactic sources and showed that the I3.3/I11.3 ratio was not varying significantly, while the 
I6.2/I11.3 ratio was varying by a factor of 5. In addition. Fig. [T71 shows that the effect of 
extinction do not account for our observed trends. Indeed, the bottom panel of Fig. [T71 
shows that the 8.6 /xm and 11.3 //m bands would be correlated if extinction was the main 
cause of variation of the aromatic bands. Consequently, our work shows that the variation 
of the PAH band ratios, throughout different environments, and at different spatial scales, is 
primarily controled by the ionization fraction of the molecules. In addition, it appears that 
the mixture of PAH molecules, in all these environments, is remarkably universal. 

We emphasize that our sample contains only a few AGNs. Therefore, our results re- 
fle ct primarily t he pr operti es of star forming re gions. Observations of AGN environments 
by ISmith et al.l (120071 ) and iKaneda et al.l (120071 ) show PAH properties that differ from our 
trends, suggesting that the AGN could alter the PAH composition. 



5.2. Relating the Band Ratios to the Physical Conditions 



In §5.11 we showed that the band ratio variations were essentially due to a variation of the 
fraction of ionized PAHs in the beam. The mid-IR feature ratios are therefore related to the 
physical conditions where the emission is originating. Ind eed, the fract ion of ionized PAHs, 
in a given region, depends on the quantity Go/ueX y^T^ ( lTielensll2005l . for a review), where 
Go, is the UV radiation field density, rig, the electron density, and Tgas, the gas temperature. 
The UV field density is the integration between A = 0.09 /xm and A = 0.2 /xm of the 
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monochromatic mean intensity, Jx^ and is normalised to the solar neighborhood value: 



47r Jx dA 



0.09 /im 



1.6 X 10-6 Wm-2 



(3) 



We emphasize that, although §5.1 showed that the intensity of the ISRF was not responsible 
for the observed variations of the aromatic band ratios, the intensity of the ISRF is likely 
to vary significantly within and among our sources. In this section, we use the most well- 
studied Galactic regions of our sample, in order to derive an empirical relation between the 
I6.2/I11.3 band ratio and the quantity Go/'^e x y/T^s- 



Region 


Go 


^ga. [K] 




cm ^ 


Reference 


NGC 2023 


1.5 X 10^ 


750 




10^ 




Steiman-Cameron et al. ("1997) 


NGC 7027 


6 X 10^ 


2000 




10^ 




Justtanont et al. ("1997) 


Orion bar (position 4) 


4 X 10^ 


500 


5 X 10^ 




Tauber et al. ^1994) 



Table 3: Physical conditions in select regions. The last columns refers to the study where 
these quantities have been derived. 



Bregman fc Temil ( 120051 ) gave such a calibration, from the spectro-images of three re- 
flection nebulae. In order to derive the physical conditions in their nebulae, they made 
several hypotheses. They (i) assumed a constant electronic density throughout their re- 
gions, (ii) they neglected dust extinction, and (Hi) they derived the relative geometry of the 
star and nebula from the profile of the scattered light. In our analysis, we have adopted a 
slightly diflFerent strategy. We have made the conscious decision to include only the sources 
for which the physical conditions have been reliably determined, using PDR models, con- 
strained by several gas lines. The sources of our sample which have been m odeled in detail are 
NGC 2023, NGC 7027 and position 4 of the Orion bar dRoche et al.lll989l ). Their properties 
are summarized in Table [31 

We adopt the values given in Table [3] and assume that the electrons result from photo- 
ionization of carbon and that all the gas phase carbon is ionized in PDRs, so that: 



C 



1.6 X 10"^ nn, 



(4) 



where we have adopted the carbon abundance bv ISofla et al.l (120041 ). For NGC 2023 and 
NGC 7027, we do not have reliable spatial information, hence we will only use the global 
values of the ratios. In the Orion bar, we study the spatial variations by proceeding as 
follows. 



-36- 



Orion Bar (Spline method) 



V 0.10 



Edge-on PDR 





10" 



H II region 
cavity 



10" 



CO 



Orion Bar (Lorentzian method) 



V 0.10 



1 



Edge-on PDR 



H II region 
cavity 



4i. 



WMW. 



% 



m - 



10" 



10" 



Icont [Wm" 



Fig. 18. — Top panels: variation of the color of the VSG continuum (ratio of the flux at 8 
and 14 //m) as a function of the intensity of the VSG continuum, inside the Orion bar. These 
figures show that the color of the VSG continuum is virtually independent of its intensity, 
since these grains are stochastically heated. Bottom panels: variation of the PAH-to-VSG 
ratio as a function of the intensity. These figures show that the relative strength of the PAHs 
in the PDR decreases slightly toward the Hll region (when Icont increases), as an effect of 
their photo-depletion. 
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Fig. 19. — Theoretical modeling of the continuum intensity as a function of the UV radiation 



field. Go, in black. I/cont/Mdust is the luminosity of the graphite component of IZubko et al. 



( I2OO4I ). integrated between 10 and 16 /xm, and normalized to the dust mass. The dotted 
line varies as Go- The dashed-dotted line is a power-law fit to the black line, in the range 
10^ < Go < lOl 
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We exclude the region north of the bar (see Fig. HOj), which is the location of the Hll 
region cavity. Indeed, the PAH-to-continuum ratio is very low there, thus the fits are 
uncertain. Moreover, the spherical geometry of the region causes confusion along the 
lines of sight. To the contrary, the region south of the bar, is an edge-on PDR, where 
the conditions are believed to be roughly homogeneous along each line of sight. Indeed, 
it is observed to be wel l stratified in the main PDR tracers (CO data as well as H2 



Cl, Cll and PAH data: iTielens et al.lll993l : iTauber et al.lll994l : IWvrowski et all 119971 : 
and HCN data bv lPuente et al.lll996l ). 



2. We assume that the column density is homogeneous in the edge-on PDR, so that 
we can use the intensity of the VSG continuum, in order to scale the intensity of 
the UV radiation field, Go- Indeed, the top panels of Fig. [18] show that the VSGs 
are stochastically heated. Therefore, their intensity scales with the radiation field 
density. Moreover, the intensity of the VSG continuum is probably a better tracer of 
the radiation density than the PAHs, since the latter are gradually depleted in direction 
of the Hll region (bottom panels of Fig. [18]). Fig. [T9l illustrates this fact by showing 
the theoretical variation of the integrated intensity of the VSG continuum between 10 
and 16 /im, as a function of Gq. To perform this simulation, we use d the radiation field 



defined in Eq. ([31), and the Galactic graphite size distribution of IZubko et al.l (120041 . 
with solar abundance constraints). For Gq < 10^, the mid-IR continuum intensity 
scales perfectly with Gq, since the grains are out of thermal equilibrium. At higher Gq^ 
the relation between Gq and /cont becomes non-linear due to temperature effects. In 

the range 10^ ^ Gq ^ 10^, which are the conditions in the Orion bar, the UV radiation 

1/13 

field is very well approximated by Go oc I^^^^ . This method has the advantage of 
taking into account the actual variations of Gq, with both the distance from the star 
cluster, and the dust extinction. 

3. In order to improve the signal-to-noise ratio of our spectra, we average the pixels into 
30 bins of energy (/cont)- The two top panels of Fig. [20] show the variation of the 
I6.2/I11.3 ratio, as a function of the continuum intensity, with the two methods. 

4. Finally, we assume a constant electron density and temperature throughout the PDR. 



The two bottom panel s of Fig. [201 shows t he fin al calibration of the band ratio, with the two 
methods. As found by lBregman fc Temil (120051 ) . the ratio I6.2/I11.3 (or I7.7/I11.3) increases in 
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the range 10^ ^ Go/ue ^ 10^. The grey stripes are the foUowing hnear fits: 

/3040 + 0.53 ± 0.10 {Spline method), 



l6.2 / Go Tu 



I11.3 \ ne [cm-3] V 103 K 



^gas 



l6.2 / Go To 



I11.3 \ He [cm-3] V 103 K 



^gas 



(5) 

/1990 + 0.26 ± 0.16 {Lorentzian method), 



G T 

valid in the range 400 < r^^\ < 4000. 

He cm~'^ V 10*^ K 



5.3. Diagnosing Galaxies Using PAH Band Ratios 

In principle, Eq. ([5j) can be used to derive the ratio Go/rie x a/T^ from any mid-IR 
spectrum. However, if several regions with different physical conditions are encompassed 
within the spectrograph's beam, the situation is more complex. The band ratio is then 
averaged over these various environments and its global value does not reflect a single Go/rieX 

In the case of the global spectrum of a star forming region, we can expect the I7.7/I11.3 
ratio to depend on the age and initial mass function of the stellar cluster, that will determine 
the Go at the edge of the PDR, as well as on the geometry of the ISM (dumpiness, densities), 
that will determine the transfer of the UV light into the PDR, and the electron density. If this 
hypothesis is correct, the value of I7.7/I11.3 can potentially be used in combination with other 
PDR tracers, like the [Clljiss^m, [Oijes/im, [Si ii]34.82/im and H2 mid-IR rotational lines, to 
determine the geometry and physical conditions of an unresolved star forming region. Fig.[2T] 
demonstrates the relation between the observed band ratios and the averaged Go/^e x \/^gas 
quantity. 

It has long been surmised that the IR emission features provide a clear mid-IR signature 



(Genzel et al. 


1998: 


Peeters et al. 


2004b) 



study has extended this by developing the observed 6.2 to 11.3 /xm band ratio as a quantita- 
tive tool to probe the physical conditions (e.g. Go/rie x ^/T^) in the emitting regions. We 
expect that this study will be of fundamental value for the interpretation of Spitzer data as 
well as future Sofia and JWST observations of galaxies in the nearby and early universe. 
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Fig. 20. — Empirical calibration of the band ratio. The two upper panels show the variation 
of the band ratio inside the Orion bar, as a function of the VSG intensity. The black error 
bars correspond to the edge-on PDR, and the grey error bars, to the H ll region cavity. The 
two lower panels show the final calibration, including NGC 2023 and NGC 7027. We kept 
only the edge-on PDR measurements, in the Orion bar. 



-41 - 



Spline method Lorentzian method 




Fig. 21. — Determination of Go/ue x ^/T^Q for the sources in our sample. These figures are 
similar to Fig. [3l where the I6.2/I11.3 has been converted into G^jn^ x ^JT^^ using Eq. ([5]). 
The same symbol conventions are adopted as in Fig. [3l The small additional symbols are 
the spatial variations within M 82. 

6. SUMMARY AND CONCLUSION 

We presented the results of a systematic study aimed at understanding the main proper- 
ties of the mid-IR features in different astrophysical environments. It is based on observations 
of Galactic regions and galaxies of various types, observed by the satellites ISO and Spitzer. 
We have developed two distinct methods of spectral decomposition with diflFerent hypothe- 
ses, in order to test the robustness of our trends, and overcome eventual biases. These two 
methods have shown similar trends between the physical quantities that we have studied. 
We explored the variations of the different mid-IR features between the wavelengths 5 and 
16 //m, among integrated spectra of galaxies, and inside galaxies and Galactic regions. Our 
main results are the foUowings. 

1. We find that the 6.2, 7.7, 8.6 /xm features are essentially tied together, while the ratio 
of these bands to the 11.3 /xm feature can vary by one order of magnitude, in our 
sample. These variations are seen both inside individual sources (like M 82, M 17, 
the Orion bar, 30 Doradus, etc.), as well as among integrated spectra. In general, the 
I6.2/I11.3 ratio is spatially correlated with the power radiated by the PAHs. It indicates 
that the ratio I6.2/I11.3 is higher in regions of intense star formation. 
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2. With the help of a stochastic heating model and realistic absorption efficiencies, we 
show that the variations of the mid-IR features are essentially due to the variation of 
the fraction of ionized PAHs. We conclude that the properties of the PAHs, throughout 
our sample, are remarkably universal. In particular, we rule out both the modification 
of the grain size distribution, and the extinction by the 9.7 /xm silicate feature, as 
an explanation of these variations. Indeed, we show that a modification of the size 
distribution could explain the observed variation of the I6.2/I11.3 ratio. However, it 
would cause the I6.2/I7.7 and I6.2/I8.6 (and I3.3/I11.3) ratios to vary significantly and 
to be correlated with the I6.2/I11.3, contradicting our observations. Similarly, a deep 
absorption by the silicate feature would vary the I6.2/I11.3 ratio, but it would decouple 
the 8.6 fim from the 6.2 and 7.7 /xm features at the same time. 

3. The universality of the properties of the PAHs and the fact that the band ratios are 
mainly sensitive to the charge of the molecules allow us to use these band ratios as a 
tracer of the physical conditions inside the emitting region. Using a few well-studied 
Galactic regions (including the spectral image of the Orion bar), we give an empirical 
relation between I6.2/I11.3 and the ratio Go/rie x ^/T^. 

4. In the case where several regions with different physical conditions are integrated, the 
band ratios are dependent on the morphology of the ISM, as well as on the evolutionary 
properties of the illuminating stellar cluster. Thus, we find that the I6.2/I11.3 band ratios 
in the star forming nuclei of M 51 and M 82 are very similar to those in the galactic 
compact Hll regions, IR 15384 and IR 18317, while the halo of M 82 and the spiral 
arms of M 51 are similar to those in the Galactic refiection nebula NGC 2023. These 
differences in band ratio refiect differences in the ionization over recombination rate 
and hence trace back to variations in the ratio of the ionizing radiation field to the 
electron density. 



A. ADDITIONAL TABLE AND FIGURES 

Notes on Individual Sources 

M 82. (Figs. [THE]) Fig. [71 shows that the correlations inside M 82 are remarkably well ar- 
ticulated. The correlation coefficients (Table [21) are 0.97, for the correlations between 
I6.2/I11.3 and I7.7/I11.3. The band ratios inside M 82 follow the general trends observed 
for the integrated spectra very well. The centroid of the 7.7 /xm feature is essentially 
fixed throughout M 82 (Table [2]). The ratio I6.2/I11.3 follows the spatial distribution 
of the PAH emission very well, which is maximum in the central bar, where the star 
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formation occurs. With good confidence (S/N > 6), we detect low ratios in the outer 
regions. 

IC 342. (Figs.[22H23D The morphology of the PAH emission (Fig.[37D follows the ^^CO(J=1^0)2.6mm 
emission which traces out the nuclear ring and ridge associated with the stellar bar 



Sakamoto et al.l (119991 ). The spectra of this galaxy are much noisier than those of M 82. 
The correlations are very well defined, with a linear correlation coefficient of 0.8 (Ta- 
ble [2]). The deviation of the relation involving the 8.6 /xm feature with the Lorentzian 
method (middle right panel of Fig. [22j), compared to the general correlation, is due to 
a systematic over-estimation of the attenuation, in the central region. The average 
centroid of the 7.7 fim feature is identical for IC 342 and M 82 ((A7.7) 7.68 //m; 
Table Ej). 

M 51. (Figs. [TOlfTTI) This galaxy is more quiescent than M 82 and IC 342. Its band ratios 
show less variation (Fig. [10]) throughout the galaxy, although they are in good agree- 
ment with the general correlations of Fig. [3l The correlation coefficients are 0.6 
(Table [2j) slightly lower than those for M 82 and IC 342, essentially because of the 
lower signal-to-noise ratio of the spectra. The spatial distribution of the I6.2/I11.3 ratio 
(Fig. [H]) follows the PAH emission remarkably well enhanced along the spiral arms 
and the bulge. This band ratio reaches its maximum on the central ring and decreases 
somewhat in the center. 

M 83. (Figs. [2511261) M 83 is a spiral galaxy which is similar in several aspects to M 51. Its 
band ratios (Fig. [25]) and the spatial distribution of its PAH emission (Fig. [38]) have 
the same properties qualitatively. The average centroid of the 7.7 /xm is identical to 
that for M 51 and M 83 ((A7.7) ^ 7.72 /xm; Table [2]). 

30 Doradus. (Figs. [2811291) The band ratio correlations measured inside 30 Doradus are sig- 
nificantly dispersed (Fig. [28]). The results of the 2 different feature extraction methods 
are not particularly consistent. This is certainly due to the very low PAH-to-continuum 
ratio. Indeed, Fig. [391 shows that the spatial distribution of the PAH emission, obtained 
with the Lorentzian method^ is different from that obtained with the Spline method. In 
particular, the right panel of Fig. [391 shows 3 bright spots along the bar of 30 Doradus, 
that are not seen on the left panel. After verification, these are fits of such low PAH- 
to-continuum ratio pixels that the widths of the PAH bands are not constrained at all. 
The Lorentzian method determines very wide bands whose integrated intensities are sig- 
nificantly overestimated. For the same reason, the highest values of the I6.2/I11.3 ratio, 
with the Spline method {c^ 2; Fig. [39]), correspond to low values, with the Lorentzian 
method {< 0.3). 
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M 17. (Figs. [3^11351) We encounter the same trouble in the Hll region part of M 17 (upper 
left corner of the images in Fig. [H]), as in 30 Doradus. For the same reasons, the 
I6.2/I11.3 ratio is maximum in this region with the Spline method {0^ 2.5; Fig. [¥TI) . 
while it is minimum with the Lorentzian method (< 0.7). However, the PDR part has 
relatively high PAH-to-continuum ratios. The band ratio correlations are well defined 
inside M 17, except for 8.6 /xm, with the Spline method (Fig. [3^ . 

The Orion bar. (Figs. [3T1I321) The PAH-to-continuum ratio is on average higher in this 
region than in M 17. Thus the fits are better, even along the ionization front. The 
correlation between band ratios (Fig. [3l]) are relatively good with linear correlation 
coefficients around 0.9 (Table [2j). The plot of the consistency between the two methods 
(Fig. [33]) shows a few outlying points. These fits are those of the pixels located inside 
the Hll region (upper part of the images in Fig. HO]). The spatial distribution of the 
I6.2/I11.3 band ratio decreases on average from the north, where the ionizing stars are 
located, to the south with the two methods. The average centroid of the 7.7 /xm feature 
is similar in M 17 and the Orion bar ((A7.7) 7.68 /xm; Table [2j). 



Table 4. Results of the spectral decomposition. 



Name Is. 2 I7.7 ^8.6 111. 3 IpAH Icont A7.7 Relia- 

[10-15 Wm-2] [10-15 Wm-2] [lO-^^Wm-^] [lO-^^Wm-^] [lO-^^Wm-^] [IQ-l^Wm-^] [^m] bility 



Haro 11 


S 


0.90 ± 0.09 


1.81 ± 0.18 


0.23 ± 0.02 


3.78 ± 2.15 


3.84 ± 2.15 


51.1 ± 5.1 






Low 




C 


1.59 ± 0.16 


5.23 ± 0.52 


0.32 ± 0.03 


1.28 ± 0.13 


11.2 ± 0.6 


59.8 ± 6.0 


7.634 


± 0.025 


Low 


SMC Bl#l 


S 


16.2 ± 2.4 


24.1 ± 3.6 


3.09 ± 0.46 


16.0 ± 2.4 


73.8 ± 35.7 


76.8 ± 11.5 






High 




C 


33.8 ± 5.1 


70.2 ± 1.1 


16.1 ± 2.4 


51.0 ± 7.6 


205 ± 14 


15.1 ± 2.3 


7.647 


± 0.025 


High 


NGC 253 


S 


148 ± 4 


293 ± 11 


28.5 ± 3.7 


105 ± 4 


604 ± 346 


(1.74 ± 0.09) X 10^ 






High 




C 


532 ± 28 


(1.54 ± 0.06) X 10^ 


36.4 ± 1.7 


313 ± 13 


(3.42 ± 0.08) X 10^ 


(1.49 ± 0.07) X 10^ 


7.707 


± 0.025 


High 


NGC 253 p 


S 


48.5 ± 2.2 


98.0 ± 6.3 


10.2 ± 0.2 


29.1 ± 1.9 


195 ± 114 


823 ± 44 






Low 




C 


15.1 ± 8.3 


481 ± 22 


113 ± 6 


114 ± 6 


(1.09 ± 0.03) X 10^ 


(1.00 ± 0.05) X 10^ 


7.702 


± 0.025 


High 


NGC 253 e 


S 


102 ± 3 


199 ± 8 


19.1 ± 0.6 


78.6 ± 2.6 


420 ± 238 


974 ± 66 






High 




C 


324 ± 22 


(1.00 ± 0.05) X 10^ 


228 ± 13 


225 ± 11 


(2.21 ± 0.06) X 10^ 


689 ± 40 


7.711 


± 0.025 


High 


SMC N 66 


S 


10.3 ± 0.5 


14.0 ± 0.6 


0.34 ± 0.05 


6.09 ± 0.11 


37.0 ± 19.1 


290 ± 41 






Low 




C 


82.7 ± 20.6 


99.7 ± 19.1 


57.5 ± 12.9 


64.5 ± 8.0 


340 ± 34 


203 ± 28 


7.619 


± 0.025 


Low 


NGC 520 


S 


8.38 ± 0.37 


19.6 ± 0.7 


1.00 ± 0.05 


4.50 ± 0.10 


36.6 ± 22.0 


45.4 ± 4.3 






High 




C 


38.9 ± 2.9 


107 ± 5 


29.2 ± 1.9 


26.3 ± 1.7 


240 ± 7 


72.3 ± 5.2 


7.748 


± 0.025 


High 


NGC 613 


S 


7.92 ± 0.26 


13.5 ± 1.0 


1.49 ± 0.01 


4.89 ± 0.05 


30.2 ± 16.5 


74.0 ± 31.9 






Low 




C 


19.3 ± 9.2 


40.4 ± 10.4 


7.41 ± 3.16 


15.9 ± 3.6 


112 ± 18 


55.4 ± 26.4 


7.693 


± 0.025 


Low 


NGC 613 p 


S 


4.03 ± 0.12 


8.35 ± 0.48 


0.72 ± 0.05 


3.88 ± 0.20 


18.3 ± 10.1 


35.9 ± 10.0 






High 




C 


12.4 ± 3.5 


32.4 ± 5.2 


6.94 ± 1.67 


12.7 ± 1.6 


79.8 ± 7.5 


33.9 ± 10.8 


7.703 


± 0.025 


High 


NGC 891 


S 


93.1 ± 9.3 


181 ± 18 


17.0 ± 1.7 


55.0 ± 5.5 


361 ± 212 


304 ± 30 






High 




C 


207 ± 21 


605 ± 61 


124 ± 12 


148 ± 15 


(1.34 ± 0.07) X 10^ 


87.1 ± 8.7 


7.705 


± 0.025 


High 


NGC 1068 


S 


55.5 ± 5.9 


126 ± 35 


20.1 ± 6.1 


59.3 ± 33.6 


287 ± 152 


(4.80 ± 0.37) X 10^ 






Low 




C 


181 ± 10 


586 ± 85 


131 ± 16 


216 ± 11 


(1.33 ± 0.09) X 10^ 


(4.75 ± 0.23) X 10^ 


7.706 


± 0.025 


Low 


NGC 1068 p 


S 


< 10.2 


< 45 


5.1 ± 4.8 


< 31 


573 ± 274 


(3.59 ± 0.34) X 10^ 






Low 




C 


209 ± 15 


352 ± 80 


86.0 ± 15.0 


97.3 ± 6.0 


902 ± 83 


(3.52 ± 0.19) X 10^ 


7.699 


± 0.025 


Low 


NGC 1068 e 


S 


48.9 ± 7.1 


101 ± 35 


14.3 ± 6.0 


51.1 ± 36.4 


237 ± 125 


(1.01 ± 0.30) X 10^ 






Low 




c 


313 ± 43 


537 ± 202 


108 ± 40 


115 ± 13 


(1.28 ± 0.21) X 10^ 


859 ± 136 


7.699 


± 0.025 


High 


NGC 1097 


s 


27.7 ± 0.1 


50.6 ± 0.01 


5.25 ± 0.15 


27.0 ± 0.8 


116 ± 64 


229 ± 32 






High 




c 


56.4 ± 6.9 


166 ± 12 


31.0 ± 3.1 


61.7 ± 4.7 


410 ± 18 


127 db 18 


7.707 


± 0.025 


High 


NGC 1097 p 


s 


18.7 ± 0.15 


37.7 ± 1.2 


4.25 ± 0.32 


19.1 ± 1.0 


82.8 ± 46.4 


156 ± 13.8 






High 




c 


42.8 ± 3.2 


134 ± 7 


29.5 ± 2.0 


48.4 ± 3.1 


319 ± 9 


88.9 ± 8.0 


7.703 


± 0.025 


High 


NGC 1140 


s 


0.92 ± 0.05 


1.31 ± 0.01 


0.07 ± 0.02 


0.62 ± 0.01 


3.34 ± 1.73 


6.96 ± 2.16 






Low 




c 


1.63 ± 0.66 


3.65 ± 0.81 


0.45 ± 0.13 


1.45 ± 0.17 


8.29 ± 1.17 


5.57 ± 1.71 


7.684 


± 0.025 


Low 


NGC 1365 


s 


38.0 ± 0.5 


76.5 ± 2.4 


9.21 ± 0.27 


34.9 ± 1.3 


177 ± 94 


367 ± 31 






High 




c 


89.3 ± 7.1 


253 ± 13 


52.7 ± 3.8 


79.5 ± 4.4 


566 ± 16 


242 ± 20 


7.710 


± 0.025 


High 


NGC 1365 p 


s 


26.7 ± 1.1 


54.1 ± 1.9 


6.22 ± 0.21 


22.3 ± 1.2 


117 ± 65 


253 ± 15 






High 




c 


57.7 ± 3.1 


182 ± 8 


42.3 ± 2.2 


53.2 ± 2.5 


148 ± 10 


71.1 ± 16.2 


7.713 


± 0.025 


High 


NGC 1365 e 


s 


11.4 ± 0.02 


22.4 ± 1.0 


2.98 ± 0.10 


12.6 ± 0.7 


60.5 ± 29.6 


113 ± 24 






High 




c 


18.3 ± 4.9 


61.2 ± 7.6 


11.7 ± 2.4 


30.8 ± 3.1 


148 ± 10 


71.1 ± 16.2 


7.708 


± 0.025 


High 


IC 342 


s 


36.8 ± 0.7 


62.1 ± 1.6 


9.69 ± 0.27 


26.7 ± 0.8 


141 ± 78 


270 ± 15 






High 




c 


62.2 ± 3.9 


182 ± 8 


48.9 ± 2.6 


60.7 ± 2.6 


441 ± 10 


158 ± 8 


7.686 


± 0.025 


High 


IC 342 p 


s 


16.5 ± 0.6 


28.3 ± 1.1 


3.90 ± 0.16 


11.3 ± 0.5 


62.9 ± 34.9 


135 ± 6 






High 




c 


30.4 ± 1.6 


78.0 ± 3.2 


15.2 ± 0.8 


26.9 ± 1.2 


188 ± 4 


87.5 ± 3.9 


7.678 


± 0.025 


High 


IC 342 e 


s 


20.2 ± 0.3 


33.8 ± 0.7 


5.79 ± 0.14 


15.4 ± 0.4 


78.4 ± 42.7 


135 ± 12 






High 




c 


32.8 ± 3.3 


105 ± 6 


34.4 ± 2.5 


33.5 ± 1.8 


257 ± 8 


70.1 ± 5.4 


7.692 


± 0.025 


High 


NGC 1569 


s 


6.82 ± 0.03 


13.2 ± 0.1 


0.50 ± 0.08 


5.38 ± 0.12 


28.1 ± 15.8 


114 ± 24 






Low 




c 


20.2 ± 9.8 


35.7 ± 10.3 


3.21 ± 1.35 


13.9 ± 2.2 


102 ± 29 


91.4 ± 14.9 


7.669 


± 0.025 


Low 


NGC 1569 e 


s 


5.43 ± 0.03 


10.7 ± 0.1 


0.33 ± 0.08 


4.02 ± 0.10 


22.5 ± 12.7 


55.8 ± 23.9 






Low 




c 


13.3 ± 9.6 


29.5 ± 12.3 


3.06 ± 1.98 


10.3 ± 2.5 


69.7 ± 38.2 


46.2 ± 15.2 


7.679 


± 0.025 


Low 


NGC 1808 


s 


86.8 ± 8.7 


147 ± 15 


15.0 ± 1.5 


59.1 ± 5.9 


321 ± 181 


497 ± 50 






High 



or 



Table 4 — Continued 



Name Iq_2 I?. 7 ^8.6 111. 3 IpAH Icont A7.7 Relia- 

[10-15 Wm-2] [10-15 Wm-2] [IQ-l^Wm-^] [IQ-l^Wm-^] [IQ-l^Wm-^] [IQ-l^Wm-^] [^m] bility 





C 




180 ± 18 


487 ± 49 


102 ± 10 


185 ± 19 


(1.24 ± 0.06) X 10^ 




274 ± 27 


7, 


,710 




0, 


,025 


High 


Orion bar D8 


S 




338 ± 3 


632 ± 7 


77.8 ± 4.4 


301 ± 12 


(1.45 ± 0.78) X 10^ 


(9, 


,80 ± 0.11) X 10^ 












Low 




C 




590 ± 13 


912 ± 5 


73.8 ± 0.8 


501 ± 12 


(2.79 ± 0.03) X 10^ 


(9, 


,25 ± 0.08) X 10^ 


7, 


,645 




0, 


,025 


Low 


Orion bar D5 


S 




532 ± 3 


(1.05 ± 0.01) X 10^ 


131 ± 2 


470 ± 3 


(2.33 ± 1.28) X 10^ 


(5, 


,84 ± 0.13) X 10^ 












High 




C 




981 ± 18 


(2.49 ± 0.02) X 10^ 


580 ± 7 


737 ± 7 


(5.97 ± 0.04) X 10^ 


(4. 


,94 ± 0.07) X 10^ 


7, 


,670 







.025 


High 


Orion bar D2 


S 




107 ± 10 


186 ± 1 


21.0 ± 0.1 


102 ± 8 


468 ± 240 


(2, 


,67 ± 0.12) X 10^ 












Low 




C 




233 ± 10 


491 ± 8 


70.1 ± 1.6 


182 ± 5 


(1.58 ± 0.03) X 10^ 


(2, 


,18 ± 0.10) X 10^ 


7, 


,680 







.025 


High 


30 Doradus 


S 




205 ± 1 


405 ± 8 


12.4 ± 0.3 


98.2 ± 3.5 


78.4 ± 46.6 


(!■ 


,39 ± 0.06) X 10^ 












Low 




C 


(1.12 ± 0.07) X 10^ 


(1.68 ± 0.07) X 10^ 


781 ± 36 


(1.04 ± 0.04) X 10^ 


(5.73 ± 0.13) X 10^ 


(!■ 


,23 ± 0.05) X 10^ 


7, 


,616 


± 





.025 


Low 


30 Doradus p 


s 




25.3 ± 1.2 


56.0 ± 0.1 


0.73 ± 0.16 


11.0 ± 0.4 


99.7 ± 62.6 


(2, 


,84 ± 0.11) X 10^ 












Low 




c 




184 ± 23 


264 ± 18 


122 ± 8 


134 ± 6 


873 ± 36 


(2, 


,61 ± 0.10) X 10"^ 


7, 


,711 


± 





.025 




30 Doradus s 


s 




164 ± 1 


315 ± 4 


11.7 ± 0.2 


81.2 ± 2.9 


624 ± 366 


(9 


,62 zb 0.40) X 10"^ 












Low 




£ 




839 ± 61 


(1.34 zb 0.07) X 10"^ 


479 zb 26 


798 zb 33 


(4.30 zb 0.12) X 10"^ 


f8 


,49 zb 0.35) X 10"^ 


7 


,624 


-|- 





.025 


High 


NGC 2023 






54.1 ± 2.3 


125 ± 5 


8.25 zb 0.11 


58.2 zb 0.5 


263 zb 149 




263 zb 55 












Hip-h 
jnign 








120 ± 9 


284 zb 8 


QQ g _[_ 3 7 


106 zb 7 


705 zb 19 




99.7 zb 14 


7 


,743 


-|- 





.025 


Hip-h 
jnign 


He 2-10 






6.74 ± 0.67 


14.7 ± 1.5 


^.U± ZIZ u.^u 


c nn -L n KQ 

O.cyU ZIZ U.Ocy 


31.1 zb 17.4 




109 zb 11 




















14.4 ± 1.4 


Of! -1 _|_ q f! 

ou . ± ziz . u 


(J.OU ZIZ \J.KJKJ 


14.2 zb 1.4 


oc n _|_ 4 q 
00 . cy ZIZ ^ . 




125 zb 13 


7 


,662 


-|- 





.025 


xiign 


M 82 






992 ± 36 


ni -1- n n7^ v in*^ 


231 zb 8 


Rn9 -1- 99 
HZ zz 


(A no -1- 9 "^zL V 1 




7zL -1- n 9Q^ V 1 

,/^ZIZU.Zb'j A ±U 












nign 






(2.6 


8 ± 0.10) X 10^ 


(^/.DZZtU.Z/j X lU 


/I on _|_ n (T/^ V 1 n3 


/o OQ _|_ n nn^ v 1 n3 
(^Z.Zo zt u.uyj X lU 


/I 70 _|_ n nq^ v 1 n4 


[p. 


,oZ ztz U.Z41:j X lU 




697 


_|_ 




095 


High 


M 






369 ± 18 


/yu It 00 


88.4 zb 3.5 




^1 c;q _|_ n nn ^ 1 n3 
(^l.oo zt u.yu X lU 


(3, 


i«_|_n 19^ V 1 n3 
,10 ztz U.lZj X lU 












High 




c 




983 ± 42 


(z.yz it u.iij X lu 


oy I ztz zo 


'7ni _L 0'7 

^ Ul ± Z ( 


( R Q c: _|_ n 1 o^ 1 n3 
(0.00 ztz U.lzj X lU 


(3, 


n9 _L n 1 1 ^ 1 n3 
,Uz ztz U. 11 j X lU 


7, 


,695 


zb 





.025 


High 


M 82 e 


s 




819 ± 30 


(1.d4 ± U.UDj X lU 


192 zb 7 


OT —I— on 
oz / ztz zU 


/ Q on —I— 1 no \/ 1 n 3 
(o.oy ztz l.yo X lU 


(5, 


—I— n o/l^ \/ in3 
,zz ztz U.z4 j X lU 












High 




c 


(1.99 ± 0.08) X 10^ 


(5.71 ± 0.21) X 10"^ 


(1.21 ± 0.04) X 10"^ 


(1.26 ± 0.05) X 10"^ 


(1.21 ± 0.02) X 10"^ 


(2, 


,82 ± 0.13) X 10"^ 


7, 


,701 


± 





.025 


High 


iNI CjU OZOD 


s 




24.7 ± 0.3 


45.5 zb 4.2 


5.29 zb 0.02 


19.4 zb 0.2 


97.4 zb 55.6 




193 zb 27 












High 




c 




65.0 ± 8.9 


183 ± 14 


48.6 ± 5.9 


57.9 ± 4.5 


411 ± 19 




181 ± 20 


7, 


,701 


± 





.025 


High 


NGC 3256 p 


s 




6.16 ± 0.40 


11.0 ± 0.3 


1.05 ± 0.02 


4.37 ± 0.25 


23.1 ± 13.4 




52.1 ± 10.2 












High 




c 




18.3 ± 3.9 


45.4 ± 5.8 


9.88 zb 2.31 


13.6 ± 1.7 


105 ± 10 




51.7 zb 7.3 


7, 


,703 


zb 





.025 


High 


NGC 3256 e 


s 




14.5 ± 1.5 


24.5 ± 0.1 


3.19 zb 0.01 


12.3 ± 0.1 


55.7 zb 32. 1 




101 ± 26 












High 




c 




38.0 ± 9.3 


104 zb 13 


29.2 zb 6.0 


37. 1 zb 4.6 


245 zb 19 




89.8 zb 18.8 


7, 


,699 


zb 





.025 


High 


Mrk 33 


s 




2.51 ± 0.25 


5.12 ± 0.51 


0.31 ± 0.03 


3.14 ± 0.31 


11.6 ± 6.5 




27.8 ± 2.8 












High 




c 




3.88 ± 0.39 


14.2 ± 0.14 


2.35 ± 0.24 


4.94 ± 0.49 


30.4 ± 1.6 




25.4 ± 2.5 


7, 


,721 


± 





.025 


High 


Arp 299 


s 




9.42 ± 0.72 


16.8 ± 1.0 


1.93 ± 0.01 


6.10 ± 0.20 


37.2 ± 20.4 




223 ± 27 












Low 




c 




28.0 ± 6.6 


43.7 ± 5.7 


4.99 ± 0.82 


13.2 ± 1.6 


147 ± 18 




413 ± 36 


7, 


,669 







.025 


Low 


UM 448 


s 




0.83 ± 0.08 


1.39 ± 0.14 


0.27 ± 0.03 


0.83 ± 0.08 


3.59 ± 1.85 




11.4 ± 1.1 












High 




c 




1.34 ± 0.13 


4.92 ± 0.49 


0.90 ± 0.09 


1.52 ± 0.15 


10.0 ± 0.6 




9.95 ± 1.0 


7, 


,670 







.025 


High 


IR 12331 


s 




38.0 ± 4.2 


69.4 ± 1.9 


7.11 ± 1.20 


23.2 ± 0.7 


160 ± 84 




869 ± 62 












Low 




c 




115 ± 7 


247 ± 6 


34.5 ± 3.2 


38.5 ± 4.1 


624 ± 20 




765 zb 28 


7, 


,696 







.025 


High 


NGC 4945 


s 




65.7 ± 1.3 


182 ± 7 


11.6 ± 0.3 


40.7 ± 0.3 


310 ± 198 




393 ± 35 












High 




c 




260 ± 16 


(1.15 ± 0.05) X 10^ 


186 ± 11 


222 ± 12 


(2.15 ± 0.05) X 10^ 




389 ± 29 


7, 


,713 







.025 


High 


NGC 4945 p 


s 




29.6 ± 0.1 


94.2 ± 3.9 


3.72 ± 0.03 


10.0 ± 0.1 


149 ± 100 




156 ± 12 












High 




c 




322 ± 15 


(1.15 ± 0.05) X 10^ 


195 ± 11 


208 ± 13 


(2.22 ± 0.05) X 10^ 




370 ± 23 


7, 


,709 







.025 


High 


NGC 4945 e 


s 




36.5 ± 0.7 


87.7 ± 3.6 


7.89 ± 0.19 


30.7 ± 0.1 


164 ± 100 




237 ± 31 












High 




c 




90.1 ± 9.6 


413 ± 23 


87.5 ± 7.7 


99.1 ± 6.4 


826 ± 28 




136 ± 16 


7, 


,713 







.025 


High 


Centaurus A 


s 




17.0 ± 0.2 


26.0 ± 0.1 


2.66 ± 0.14 


13.6 ± 0.2 


61.1 ± 34.0 




219 ± 33 












Low 




c 




30.2 ± 5.8 


112 ± 13 


27.2 ± 1.3 


50.3 ± 5.0 


282 ± 16 




219 ± 30 


7, 


,696 


± 





.025 


High 


Centaurus A e 


s 




14.3 ± 0.1 


22.7 ± 0.1 


2.07 ± 0.12 


10.7 ± 0.1 


51.3 ± 29.0 




97.2 ± 32.0 












Low 




c 




23.8 ± 7.4 


87.3 ± 14.3 


21.8 ± 5.8 


36.7 ± 5.7 


219 ± 19 




64.3 ± 20.4 


7, 


,695 


± 





.025 


Low 


M 51 


s 




64.1 ± 0.9 


123 ± 3 


12.9 ± 0.4 


57.2 ± 1.7 


264 ± 151 




435 ± 49 












High 




c 




125 ± 10 


398 ± 20 


76.0 ± 5.1 


131 ± 8 


843 ± 25 




251 ± 35 


7, 


,716 







.025 


High 



Table 4 — Continued 



Name 




le 2 


I7 7 


Is 6 


111 3 


IpAH 


Icont 


A7 7 




Relia- 






[10~^^ W m~2] 


[10~^^ Wm~2] 


[10~^^ Wm~2] 


[IQ-l^ Wm-2] 


[IQ-l^ Wm-2] 


[IQ-l^ Wm-2] 


[uml 




bility 


M 51 p 


S 


22.6 ± 0.5 


42.7 ± 0.8 


4.93 ± 0.1 


21.0 ± 0.5 


93.6 ± 52.9 


161 ± 17 






High 




C 


46.0 ± 3.9 


145 ± 8 


28.8 ± 2.1 


49.5 ± 2.7 


326 ± 10 


82.4 ± 10.6 


7.721 ± 0, 


,025 


High 


M 51 e 


S 


34.9 ± 0.4 


68.4 ± 1.0 


6.76 ± 0.19 


30.7 ± 0.8 


145 ± 83 


231 ± 41 






High 




C 


61.6 ± 8.4 


212 ± 15 


41.7 ± 4.4 


71.5 ± 5.4 


436 ± 19 


141 ± 32 


7.713 ± 0, 


,025 


High 


M 83 


S 


152 ± 3 


291 ± 9 


33.9 ± 1.4 


130 ± 5 


656 ± 357 


901 ± 174 






High 




C 


304 ± 22 


810 ± 48 


147 ± 14 


343 ± 30 


(1.79 ± 0.06) X 10^ 


622 ± 154 


7.708 ± 0, 


,025 


High 


M 83 p 


S 


60.4 ± 2.8 


117 ± 5 


14.2 ± 0.5 


52.9 ± 2.7 


260 ± 143 


458 ± 28 






High 




C 


163 ± 8 


430 ± 18 


103 ± 5 


119 ± 5 


931 ± 22 


311 ± 19 


7.707 ± 0, 


,025 


High 


M 83 e 


S 


99.4 ± 1.9 


187 ± 6 


21.1 ± 0.9 


84.3 ± 3.2 


426 ± 230 


501 ± 156 






High 




C 


163 ± 16 


540 ± 40 


105 ± 14 


247 ± 28 


(1.21 ± 0.06) X 10^ 


180 ± 94 


7.712 ± 0, 


,025 


High 


Circinus 


S 


79.1 ± 0.5 


211 ± 6 


24.8 ± 0.6 


89.4 ± 5.5 


421 ± 244 


(22.8 ± 1.0) X 10^ 






Low- 




C 


347 ± 18 


(1.05 ± 0.04) X 10^ 


173 ± 8 


494 ± 20 


(2.55 ± 0.05) X 10^ 


(3.58 ± 0.14) X 10^ 


7.684 ± 0, 


,025 


High 


Circinus e 


S 


64.2 ± 0.4 


161 ± 4 


20.4 ± 0.4 


70.5 ± 5.1 


330 ± 188 


945 ± 67 






Low 




C 


261 ± 19 


774 ± 38 


187 ± 11 


385 ± 19 


(1.95 ± 0.05) X 10^ 


(1.13 ± 0.06) X 10^ 


7.688 ± 0, 


,025 


High 


Arp 220 


S 


0.48 ± 0.03 


6.33 ± 0.37 


(3.63 ± 1.72) X 10-2 


1.50 ± 0.02 


11.5 ± 7.1 


531 ± 214 






Low 




C 


< 59 


68.5 ± 23.7 


43.7 ± 32.4 


36.2 ± 21.1 


259 ± 76 


258 ± 62 


7.704 ± 0, 


,025 


Low 


Arp 220 p 


S 


0.63 ± 0.14 


4.00 ± 0.05 


< 4.4 X 10-2 


0.90 ± 0.01 


6.64 ± 4.23 


34.9 ± 9.9 






Low 




C 


31.6 ± 21.7 


68.7 ± 17.6 


78.3 ± 38.7 


57.0 ± 28.8 


334 ± 73 


247 ± 41 


7.715 ± 0, 


,025 


Low 


IR 15384 


S 


315 ± 2 


608 ± 6 


51.4 ± 2.8 


142 ± 3 


(1.22 ± 0.70) X 10^ 


(3.60 ± 0.07) X 10^ 






High 




C 


711 ± 13 


(2.13 ± 0.02) X 10^ 


304 ± 0.07 


305 ± 4 


(4.32 ± 0.03) X 10^ 


(2.89 ± 0.03) X 10^ 


7.674 ± 0, 


,025 


High 


NGC 6240 


S 


4.22 ± 0.06 


9.05 ± 0.01 


0.91 ± 0.05 


3.33 ± 0.04 


18.2 ± 10.6 


51.2 ± 6.1 






High 




c 


14.1 ± 2.6 


43.2 ± 4.0 


6.42 ± 1.05 


17.6 ± 1.5 


104 ± 6 


96.0 ± 8.5 


7.725 ± 0, 


,025 


High 


NGC 6240 p 


s 


2.88 ± 0.13 


5.78 ± 0.17 


0.49 ± 0.05 


1.94 ± 0.02 


11.4 ± 6.8 


32.6 ± 3.0 






High 




c 


10.3 ± 1.3 


30.1 ± 2.0 


3.78 ± 0.43 


9.21 ± 0.63 


66.3 ± 2.8 


59.4 ± 4.4 


7.722 ± 0, 


,025 


High 


IR 18317 


s 


292 ± 17 


527 ± 4 


44.7 ± 0.1 


122 ± 5 


(1.07 ± 0.62) X 10^ 


(3.70 ± 0.87) X 10^ 






Low 




L 


671 ± 14 


(1.99 ± 0.02) X 10^ 


256 ± 3 


275 ± 4 


(4.18 ± 0.03) X 10^ 


(2.93 ± 0.05) X 10^ 


7.681 ± 0, 


,025 


High 


NGC 6946 


S 


81.4 ± 2.8 


156 ± 6 


16.0 ± 0.5 


62.1 ± 2.4 


322 ± 188 


433 ± 47 






High 




C 


173 ± 9 


527 ± 24 


88.7 ± 5.6 


165 ± 11 


(1.10 ± 0.03) X 10^ 


192 ± 24 


7.713 ± 0, 


,025 


High 


NGC 6946 p 


S 


24.8 ± 1.8 


51.7 ± 3.3 


4.62 ± 0.04 


17.6 ± 1.3 


100 ± 60 


162 ± 10 






High 




C 


77.9 ± 4.1 


211 it 8 


42.3 ± 2.1 


58.4 ± 4.3 


462 ± 11 


85.5 ± 4.6 


7.714 ± 0, 


,025 


High 


NGC 6946 e 


s 


56.8 ± 1.8 


105 ± 4 


11.4 ± 0.4 


44.6 ± 1.8 


222 ± 128 


273 ± 39 






High 




c 


105 ± 6 


327 ± 17 


52.6 ± 4.0 


119 ± 9 


687 ± 21 


122 ± 22 


7.712 ± 0, 


,025 


High 


NGC 7027 


s 


1.00 ± 0.10 


1.82 ± 0.18 


0.14 ± 0.02 


1.52 ± 0.15 


6.84 ± 0.68 


33.4 ± 3.3 






Low 




c 


1.38 ± 0.14 


4.05 ± 0.41 


0.48 ± 0.05 


4.34 ± 0.43 


17.0 ± 0.7 


27.0 ± 2.7 


7.774 ± 0, 


,025 


High 


IR 22308 


s 


102 ± 2 


236 ± 3 


25.0 ± 1.0 


85.1 ± 1.3 


485 ± 273 


906 ± 75 






High 




c 


191 ± 5 


572 ± 10.9 


100 ± 2 


153 ± 10 


(1.31 ± 0.02) X 10^ 


629 ± 26 


7.684 ± 0, 


,025 


High 


IR 23030 


s 


75.4 ± 8.7 


136 ± 3 


17.4 ± 0.4 


54.0 ± 2.3 


301 ± 166 


(1.36 ± 0.06) X 10^ 






Low 




c 


135 ± 5 


188 ± 5 


21.5 ± 0.6 


129 ± 5 


736 ± 15 


(1.13 ± 0.03) X 10^ 


7.649 ± 0, 


,025 


High 


IR 23133 


s 


233 ± 6 


496 ± 2 


66.8 ± 0.4 


158 ± 1 


997 ± 575 


(2.88 ± 0.10) X 10^ 






High 




c 


493 ± 16 


(1.18 ± 0.02) X 10^ 


216 ± 2 


(2.71 ± 0.03) X 10^ 


(2.71 ± 0.03) X 10^ 


(2.46 ± 0.04) X 10^ 


7.674 ± 0, 


,025 


High 


IR 23128 


s 


0.82 ± 0.18 


2.43 ± 0.12 


0.27 ± 0.09 


1.14 ± 0.02 


6.35 ± 3.09 


26.4 ± 5.5 






Low 




c 


5.24 ± 1.95 


13.1 ± 2.6 


2.24 ± 0.63 


3.38 ± 0.54 


35.7 ± 5.2 


35.4 ± 5.5 


7.681 ± 0, 


,025 


Low 


IR 23128 p 


s 


0.79 ± 0.02 


1.69 ± 0.07 


0.25 ± 0.01 


0.74 ± 0.01 


4.23 ± 2.10 


11.9 ± 2.6 






Low 




c 


1.91 ± 0.56 


7.23 ± 1.06 


2.57 ± 0.56 


1.47 ± 0.20 


20.8 ± 1.6 


6.66 ± 0.90 


7.688 ± 0, 


,025 


Low 


NGC 7714 


s 


3.17 ± 0.32 


6.00 ± 0.60 


0.90 ± 0.09 


2.78 ± 0.28 


13.8 ± 7.4 


43.5 ± 4.4 






High 




c 


5.15 ± 0.51 


14.9 ± 1.5 


2.75 ± 0.28 


5.34 ± 0.53 


34.8 ± 1.7 


36.4 ± 3.6 


7.681 ± 0, 


,025 


High 



Note. — We consider that a given fit has a "high" reliability, if S/N > 6 at A = 7.7 fim, and IpAH/Icont ^ 0.3 with the Spline method, and IpAH/Icont ^ 0.5 with the Lorentzian 
method. Otherwise, the reliability is "low". 
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Fig. 22. — PAH band ratios within IC 342. The same symbol conventions are adopted as in 

Fig.m 
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Fig. 24. — Comparison between the two methods in IC 342. The same symbol conventions 
are adopted as in Fig. [71 
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M 83 (Lorentzian method) 




M 83 (Spline method) 




Fig. 25. — PAH band ratios within M 83. The same symbol conventions are adopted as in 

Fig.m 
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M 83 (Lorentzian method) M 83 (Spline method) 




Fig. 26. — PAH band ratios within M 83 (continued). 



M 83 
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Fig. 27. — Comparison between the two methods in M 83. The same symbol conventions 
are adopted as in Fig. [71 
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30 Dor (Lorentzian method) 




Fig. 28. — PAH band ratios within 30 Doradus. The same symbol conventions are adopted 
as in Fig. [3 
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30 Dor (Lorentzian method) 30 Dor (Spline method) 




Fig. 29. — PAH band ratios within 30 Doradus (continued). 
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Fig. 30. — Comparison between the two methods in 30 Doradus. The same symbol conven- 
tions are adopted as in Fig. [71 
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Orion Bar (Lorentzian method) 




Fig. 31. — PAH band ratios within the Orion bar. The same symbol conventions are adopted 
as in Fig. [3 
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Orion Bar (Lorentzian method) Orion Bar (Spline method) 




Fig. 32. — PAH band ratios within Orion bar (continued). 
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Fig. 33. — Comparison between the two methods in Orion bar. The same symbol conventions 
are adopted as in Fig. [71 
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M 17 (Lorentzian method) M 17 (Spline method) 




Fig. 35. — PAH band ratios within M 17 (continued). 




Fig. 36. — Comparison between the two methods in M 17. The same symbol conventions 
are adopted as in Fig. [71 
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IC 342: Ip,Htot &^ IpAh(6.2)/Ipah(11-3) Lorentzian 

30" 





1 1 




o 








^^^^^^^^^^^^ 




1 



50'' 45'' 
^'JSOOO 



Fig. 37. — Spatial distribution of the PAHs in IC 342. For each method, the image is the 
total PAH intensity, and the contours are the I6.2/I11.3 ratio. The shaded circle indicates the 
beam size. 
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Fig. 38. — Spatial distribution of the PAHs in M 83. For each method, the image is the 
total PAH intensity, and the contours are the I6.2/I11.3 ratio. The shaded circle indicates the 
beam size. 




Fig. 39. — Spatial distribution of the PAHs in 30 Doradus. For each method, the image is 
the total PAH intensity, and the contours are the I6.2/I11.3 ratio. The shaded circle indicates 
the beam size. 
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Fig. 40. — Spatial distribution of the PAHs in Orion bar. For each method, the image is the 
total PAH intensity, and the contours are the I6.2/I11.3 ratio. The shaded circle indicates the 
beam size. 
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Fig. 41. — Spatial distribution of the PAHs in M 17. For each method, the image is the 
total PAH intensity, and the contours are the I6.2/I11.3 ratio. The shaded circle indicates the 
beam size. 
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